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We develop a model approach to describe charge fluctuations and different charge phases in 
strongly correlated 3d oxides. As a generic model system one considers that of centers each with 
three possible valence states M '* described in frames of S = 1 pseudo-spin (isospin) formalism 
by an effective anisotropic non-Heisenberg Hamiltonian which includes two types of single particle 
correlated hopping and the two-particle hopping. We show that the coherent states provide the 
optimal way both to a correct mean-field approximation and respective continuous models to describe 
the pseudo-spin system. Simple uniform mean-field phases include an insulating monovalent M°- 
phase, mixed- valence binary (disproportionated) Af ± -phase, and mixed- valence ternary ("under- 
disproportionated" ) M 0,± -phase. We consider two first phases in more details focusing on the 
problem of electron/hole states and different types of excitons in M°-phase and formation of electron- 
hole Bose liquid in Af ± -phase. Pseudo-spin formalism provides a useful framework for revealing 
and describing different topological charge fluctuations, in particular, like domain walls or bubble 
domains in antiferromagnets. Electron-lattice polarization effects are shown to be crucial for the 
stabilization of either phase. All the insulating systems such as M°-phase may be subdivided to two 
classes: stable and unstable ones with regard to the formation of self-trapped charge transfer (CT) 
excitons. The latter systems appear to be unstable with regard to the formation of CT exciton 
clusters, or droplets of the electron-hole Bose liquid. The model approach suggested is believed 
to be applied to describe a physics of strongly correlated oxides such as cuprates, manganites, 
bismuthates, and other systems with charge transfer excitonic instability and/or mixed valence. 
We shortly discuss an unconventional scenario of the essential physics of cuprates and manganites 
that implies their instability with regard to the self-trapping of charge transfer excitons and the 
formation of electron-hole Bose liquid. 



I. INTRODUCTION 

The discovery of the high-T c superconductivity in 
doped cuprates, * observation of many unconventional 
properties in doped manganites with their colossal mag- 
netoresistance, bismuthates with high-T c 's, nickellates 
and many other oxides^ shows that we deal with a man- 
ifestation of novel strongly correlated states with a lo- 
cal charge instability, mixed valence, "metal-dielectric" 
duality, strong coupling of different (charge, spin, or- 
bital, structural) degrees of freedom and non-Landau be- 
haviour of quasiparticles. All this has generated a flurry 
of ideas, models and scenarios of the puzzling transport 
phenomena and stimulated the intensive studies of vari- 
ous correlation effects and charge transfer (CT) phenom- 
ena in strongly correlated systems derived in cither way 
from insulators unstable with regard to the CT fluctua- 
tions. Conventional approach to hotly debated strongly 
correlated 3d oxides such as cuprates, manganites, and 
many other similar systems implies making use of a Hub- 
bard model with famous Hamiltonian 

H = - ^2 *(*i)cJ<rCj> + U ^2 n ir ,n irj , (1) 

<i,j>,a i.aa' 

with competing contributions of kinetic and potential 
terms. Here c ia /cj a are creation/annihilation operators 
for low-lying antibonding 3d-0 2p hybridized orbitals. 
The two-center charge transfer integral t(ij) is often as- 
sociated with d-d transfer. Mott-Hubbard insulator is 



believed to arise from a potentially metallic half-filled 
band as a result of the Coulomb blockade of electron 
tunnelling (U ^> t) to neighboring sites£ 

Despite intense effort, the behavior of strongly corre- 
lated 3d oxides remain poorly understood and we are 
still far from a comprehensive understanding of the un- 
derlying physics. Moreover, it seems that there are miss- 
ing qualitative aspects of the problem beyond the simple 
Hubbard scenario that so far escaped the identification 
and the recognition. Firstly it concerns strong electron- 
lattice polarization effects which may be subdivided into 
electron-lattice interaction itself^*^ and a contribution 
of an electronic background that is electronic subsystem 
which is not incorporated into effective Hubbard model 
HamiltonianAi These effects are of great importance for 
the ground state electronic and crystalline structure, and 
can seriously modify the doping response of 3d oxide up 
to the crucial change of the seemingly natural ground 
state. This question has not received the attention it de- 
serves. It should be emphasized that traditional Frohlich 
approach to the electron-lattice coupling implies the de- 
scription of linear effects whereas the charge fluctuations 
in the insulator do imply strongly nonlinear electron- 
lattice coupling with the predominance of polarization 
and relaxation effects, and another energy scale. 

Electron-lattice effects may be directly incorporated 
into effective Hubbard model. Assuming the coupling 
with the local displacement (configuration) coordinate Q 
in the effective potential energy we arrive at a general- 
ized Peierls-Hubbard models From the other hand, the 
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taking account of similar effects in the kinetic energy re- 
sults in a generalized Su-Schrieffer-Heeger (SSH) models 
The correlation effect of an electronic background was 
shownSii to be of primary importance for atomic systems 
with filled or almost filled electron shells. Namely such a 
situation is realized in oxides with 2_ (2p 6 ) oxygen ions. 
In particular, the effect results in a correlated character 
of a charge transfer that seems to be one of the main 
features for 3d oxides. 

Many strongly correlated 3d oxides reveal anomalous 
sensitivity to a small nonisovalent substitution. For ex- 
ample, only 2% Sr 2+ substituted for La 3+ in La2Cu04 
result in a dramatical suppression of long-range copper 
antiferromagnetism, while it is suppressed with isovalent 
Cu 2+ substitution by Zn 2+ at a much higher concentra- 
tion close to the site dilution percolating threshold. Si- 
multaneously, the transport properties of La2- a; Sr 3 ;Cu04 
system reveal unconventional insulator-metal duality 
starting from very low dopant levels Most likely, all 
this points to a charge phase instability intrinsic for par- 
ent 214 system which somehow evolves with nonisova- 
lent substitution due to a well developed charge potential 
inhomogeneity and/or hole doping effect. The problem 
seems to be closely related with the hidden multi stability 
intrinsic to each solid&ii If the ground state of a solid 
is pseudo-degenerate, being composed of true and false 
ground states with each structural and electronic orders 
different from others, one might call it multi-stable. Be- 
low we focus ourselves on a charge degree of freedom and 
charge (in)stability, rather than orbital or spin degrees 
of freedom. As an illuminating example of such a ma- 
terial with conceptually simple but actually false ground 
state Toyozawa 11 - suggests to address the Wolfram's red, 
a quasi-one dimensional material of which the skeleton 
chain consists of alternate array: (CT~ - Pt 3+ -) 2n with 
simple (and seemingly metallic), but a false ground state. 
The real ground state is an insulator with a complicated 
structure of doubled period: (CI" - Pt 4 + - CI" - Pt 2+ -)™, 
which can be reached from the former through the Peierls 
transition with the charge density wave of large ampli- 
tude, or disproportionation like reaction. This transition 
can be considered as the condensation of self-decomposed 
self-trapped excitons spontaneously generated on all unit 
cells. 

In this connection it is worth noting the text-book ex- 
ample of BaBiOa system where we unexpectedly deal 
with the disproportionated Ba 3+ + Ba 5+ ground state 
instead of the conventional lattice of Ba 4+ cations^ 
The bismuthate situation can be viewed also as a re- 
sult of a condensation of CT excitons, in other words, 
the spontaneous generation of self-trapped CT excitons 
in the ground state with a proper transformation of lat- 
tice parameters. At present, a CT instability with re- 
gard to disproportionation is believed to be a rather typ- 
ical property for a number of perovskite 3d transition- 
metal oxides such as SrFe03, LaCu03, RNiO^, more- 
over, in solid state chemistry one consider tens of dispro- 
portionated systems^ New principles must be developed 



to treat such charge or CT unstable systems with their 
dramatical non-Fermi-liquid behavior. In particular, we 
have to change the current paradigm of the metal-to- 
insulator (MI) transition to that of an insulator-to-metal 
(IM) phase transition. These two approaches imply es- 
sentially different starting points: the former starts from 
a rather simple metallic-like scenario with inclusion of 
correlation effects, while the latter does from strongly 
correlated atomic-like scenario with the inclusion of a 
charge transfer. Electron-lattice polarization effects ac- 
companying the charge transfer appear to be of primary 
importance to stabilize either phase state. One should 
emphasize that the theoretical description of such sys- 
tems is one of the challenging problems in solid state 
physics. 

Hereafter, we develop a model approach to describe dif- 
ferent charge fluctuations and charge phases in strongly 
correlated 3d oxides with main focus on the correlated 
CT effects. As an illustrative model system we address 
a simple mixed-valence system with three possible sta- 
ble nondegenerate valent states of a cation-anionic clus- 
ter, hereafter M: M°,M ± , forming the charge (isospin) 
triplet. The M° valent state is associated with the con- 
ceptually simple one like CuO^ - in insulating copper ox- 
ides (CuO, La 2 Cu0 4 , YBa 2 Cu 3 6 , Sr 2 Cu0 2 Cl 2 ,...) or 
MnOg in manganite LaMnC>3 or BiOg - in bismuthates. 
It is worth noting that such a model is a most relevant to 
describe different cuprates where novel concepts should 
compete with a traditional Hubbard model approach in a 
hole representation implying the vacuum state formed by 
M~ (CuO^ - ) centers, and some concentration of holes. 
That is why overall the paper we refer the insulating 
cuprates to illustrate the main concepts of the approach 
developed. Our mathematics is based on the 5=1 
pseudo-spin formalism (see e.g. review article Ref lTol) 
to be the effective tool for the description of the essential 
physics both of insulators unstable with regard to the CT 
fluctuations and related mixed-valence systems. Such an 
approach provides the universal framework for a unified 
description of these systems as possible phase states of a 
certain parent multi-stable system. In addition, we may 
make use of powerful methods developed in the physics 
of spin systems. The model system of M 0,± centers is de- 
scribed in frames of S = 1 pseudo-spin formalism by an 
effective anisotropic non-Heisenberg Hamiltonian which 
includes two types of correlated two-center hopping: 

M° + AI° <-> M ± + M T andM* + M° M° + M ± , 

respectively. It should be noted that we neglect all the 
intra-center transition, including anion-cation O 2p-3d 
charge transfer. 

Our main goal is to describe different charge phases 
of the model system and a scenario of evolution of visi- 
bly typical insulator to unconventional electron-hole Bose 
liquid which reveals many unexpected properties, includ- 
ing superconductivity. The paper is organized as fol- 
lows: In Sec. II we address a metal-oxide cluster model, 
different mechanisms of correlation effects, and the ef- 



3 



fects of electron-lattice polarization. In Sec. Ill we in- 
troduce the S — 1 pseudo-spin formalism to describe 
the model mixed-valence systems. The effective pseudo- 
spin Hamiltonian and possible mean-field phase states 
of the mixed-valence systems are discussed in Sec. IV. In 
Sec.V we analyse an eh-representation of different excita- 
tions in a monovalent M° phase, discuss a CT instability, 
and nucleation of electroh-hole (EH) droplets. Electron- 
hole Bose liquid is discussed in Sec. VI. Some topological 
skyrmion-like charge fluctuations in the model mixed va- 
lence system are described in Sec. VII. Implications for 
cuprates and manganites are discussed in Sec. VIII. 



II. METAL-OXIDE CLUSTERS AND 
CORRELATION EFFECTS 
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The electronic states in strongly correlated 3d oxides 
manifest both significant correlations and dispersional 
features. The dilemma posed by such a combination is 
the overwhelming number of configurations which must 
be considered in treating strong correlations in a truly 
bulk system. One strategy to deal with this dilemma is 
to restrict oneself to small 3d-metal-oxygen clusters, cre- 
ating model Hamiltonians whose spectra may reasonably 
well represent the energy and dispersion of the impor- 
tant excitations of the full problem. Indeed, such clusters 
as CUO4 in quasi-2D cuprates, Mn06 in manganite per- 
ovskites are basic elements of crystalline and electronic 
structure. Despite a number of principal shortcomings, 
including the boundary conditions, the breaking of lo- 
cal symmetry of boundary atoms, sharing of common 
anions for nn clusters etc., the embedded molecular clus- 
ter method provides both, a clear physical picture of the 
complex electronic structure and the energy spectrum, 
as well as the possibility of quantitative modelling. Es- 
kes et ali&, as well as Ghijsen et alH have shown that 
in a certain sense the cluster calculations might provide 
a better description of the overall electronic structure of 
insulating 3d oxides than band-structure calculations. In 
particular, they allow to take better into account differ- 
ent correlation effects. 

Below, in the Section we discuss some aspects of elec- 
tronic structure, energy spectrum, and correlation effects 
for an illustrative example of Cu04 clusters embedded 
into an insulating cuprate. 



A. Electronic structure of copper-oxygen clusters 

Beginning from 5 Cu 3d and 12 O 2p atomic orbitals 
for CUO4 cluster with D±h symmetry, it is easy to form 
17 symmetrized a\ g , a2 g , b\ g , b^g, e g (gerade=even) and 
o-2u, &2u, e u (tr), e u (ir) (ungerade=odd) orbitals. The even 
Cu 3d aig(3d z 2),big(3d x 2_ y 2),b2g{3d X y),e g (3d X z,3dy Z ) 
orbitals hybridize, due to strong Cu 3d-0 2p covalency, 
with even O 2p-orbitals of the same symmetry, thus 
forming appropriate bonding j b and antibonding 7" 



FIG. 1: Model single-hole energy spectra for a C11O4 pla- 
quette with parameters relevant for a number of insulating 
cuprates. 

states. Among the odd orbitals only e u (a) and e„(7r) 
hybridize due to nearest neighbor pp overlap and trans- 
fer thus forming appropriate bonding e h u and antibonding 
e° purely oxygen states. The purely oxygen a,2 g , 02^, &2u 
orbitals are nonbonding. All "planar" O 2p orbitals in 
accordance with the orientation of lobes could be clas- 
sified as a (ai g ,bi g ,e u (cr)) or ir (d2 g , b2 g , e M (7r)) orbitals, 
respectively. 

Bonding and antibonding molecular orbitals in hole 
representation can be presented, for example, as follows 

\b\ g ) = cosa blg \b lg {3d)} +sma blg \b lg (2p)), 

K g ) = sina blg \bi g (3d)) - co$a blg \bi g (2p)), (2) 

where &i g (3d) = 3d x 2_ y 2 and \bi g (2p)} is a superposition 
of four O 2p orbitals with b\ g symmetry. 

Fig. n presents a single-hole energy spectrum for a 
CuC>4 plaquette embedded into an insulating cuprate 
like Sr2Cu02Cl2 calculated with a reasonable set of 
parameters^ For illustration we show also a step-by- 
step formation of the cluster energy levels from the bare 
Cu 3d and O 2p levels with the successive inclusion of 
crystalline field (CF) effects, O 2p-0 2p, and Cu 3d-0 
2p covalency. It is worth noting that strong Cu 3d-0 
2p overlap and covalency result in a dramatic difference 
between CuO^ - center and Cu 2+ cation as its naive ana- 
logue. 

1. Two-hole configurations for the CuO^~ center 

Starting from CuO^ - center as a realization of M° cen- 
ter we arrive at CuO|~ and CuOj - centers as realizations 
of hole M + and electron M~ centers, respectively. The 
electron CuO^" center has a filled set of Cu 3d and O 
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2p orbitals, and is often addressed to be a generalization 
of Cu 1+ ion. The hole CuO^ - center has actually two- 
hole configuration with the lowest Zhang- Rice (ZR) spin- 
singlet 1 j4i s stated formed by the interaction of three 
"covalent" configurations: (bi g ), {b\ g ), and {b\ g b\ g ), re- 
spectively, or three "purely ionic" two- hole configurations 
|d 2 ), \pd), and \p 2 ). Here, \d) = \d x i_ y i) and \p) = \p blg ) 
are the non-hybridized Cu 3d x 2_ y 2 and O 2p a orbitals, 
respectively, with bare energies e d and e p . 

Let us present a simple example of the calculation of 
the two-hole spectrum in the ZR-singlet sector. The 
matrix of the full effective Hamiltonian within the bare 
purely ionic basis set has a rather simple form 

(2e d + U d t \ 

H=\ t e d + e p + V pd t , (3) 

V o t 2e p + u;J 

where the effective Coulomb parameter for purely oxygen 
configuration incorporates both the intra-atomic param- 
eter U p and the oxygen-oxygen coupling to the first and 
second nearest neighbors, respectively 




and the following condition holds: 

U d > U p > V pd 

For reasonable values of parameters (in eV): U d = 8.5, 
U p = V pd = 1.2, e d = 0, e p = 3.0, t = t pd = 1-3 (see 
Ref. llSft we obtain for the ZR-singlet energy Ezr = 3.6, 
and its wave function 

|$i 2) ) = \b\ g ,pd) = 0.25|d 2 ) - 0.95|pd) + 0.19|p 2 ) , (4) 

where \pd) = -^(|p(l)d(2)) + |p(2)d(l))). ft reflects the 
well-known result that the ZR-singlet represents a two- 
hole configuration with one predominantly Cu 3d and 
one predominantly O 2p hole. It is worth noting that the 
hole CuO\~ center sometimes one naively associate with 
Cu 3+ ion, however, such a conclusion is a rather far from 
reality. Indeed, this center is the mixed valence one, as 
the Cu valence resonates between +2 and +3. 

The two excited states with energies Ezr + 5.2 and 
Ezr. + 6.7 eV are described by the wave functions 

|$f ) = \b 2 lgl dd) =-0.95M 2 )-0.2lH) + 0.22|p 2 ), (5) 

|<4 2) ) = \b\ gl pp) = 0.17|d 2 ) + Q.2A\pd) + 0.96|p 2 >, (6) 

respectively. Given the ZR-singlet energy one may cal- 
culate the minimal energy for M° + M° — > + 
charge transfer: 

Act = Ezr - 2E blg = (3.6 + 0.5)eV = 4.1eV, 

where the stabilization energy for the bonding b\ g state is 
simply calculated from matrix 10) at U d = U* = V pd = 0. 



It should be emphasized that this quantity plays a par- 
ticular role as the minimal charge transfer energy which 
specifies the charge transfer gap. In the general case it is 
defined as follows: 

Act — En+i + Eu-i — 2-Eat, 

or as the energy required to remove a hole from one region 
of the crystal and add it to another region beyond the 
range of excitonic correlations. The exact diagonalization 
studies for a series of clusters with different siz o 20 i 21 show 
that Act strongly diminishes with cluster size from « 
4eV for small clusters to w 2.5eV as extrapolated value 
for large clusters. 

Our simple calculation points to a significant role of 
correlation effect. Indeed, the inter-configurational cou- 
pling due to Coulomb repulsion results in a visible devia- 
tion of the two-hole ground state wave function from the 
predictions of simple model of noninteracting particles 

|(^ g ) 2 > = (-0.8|d(l))+0.6lp(l)))(-0.8|d(2)) + 0.6b(2)» 

= 0.64|d 2 ) - 0.68\pd) + 0.36|p 2 ) . (7) 

In particular, it could give rise to a strong renormaliza- 
tion of the hole transfer integrals. 

B. Electron-lattice polarization effects 

1. Correlation effects of electronic background 

The correlation problem becomes of primary impor- 
tance for atoms/ions near Coulomb instability when the 
one-electron gluing cannot get over the destructive effect 
of the electron-electron repulsion. Such a situation seems 
to realize in oxides where Hirsch et al. 6 have proposed 
an instability of 2_ (2p 6 ) electronic background. The 
main suggestion in their theory of "anionic metal" con- 
cerns the occurrence of the non-rigid degenerate struc- 
ture for a closed electron shell such as 2 ~(2p 6 ) with 
the internal purely correlation degrees of freedom. In 
other words, one should expect sizeable correlation ef- 
fects not only from unfilled 3d- or oxygen 2p shells, but 
from completely filled O 2p 6 shell! In order to relevantly 
describe such a non-rigid atomic background and its cou- 
pling to the valent hole one might use a concept of the 
well-known "shell-droplet" model for nuclei after Bohr 
and Mottelson^ In accordance with the model a set of 
completely filled electron shells which form an atomic 
background or vacuum state for a hole representation is 
described by certain internal collective degrees of free- 
dom and a number of physical quantities such as elec- 
tric quadrupole and magnetic moments. Valent hole(s) 
moves around this non-rigid background with strong in- 
teraction inbetween. Such an approach strongly differs 
from the textbook one that implies a rigid atomic orbital 
basis irrespective of varying filling number and external 
potential. 
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None of the effective many-body Hamiltonians that 
are most widely used to study the effect of electron cor- 
relation in solids such as the Hubbard model, the An- 
derson impurity and lattice models, the Kondo model, 
contain this very basic and fundamental aspect of elec- 
tron correlation that follows from the atomic analysis. 23 
The Hubbard on-site repulsion U between opposite spin 
electrons on the same atomic orbital is widely regarded 
to be the only important source of electron correlation 
in solids. It is a clear oversimplification, and we need 
in a more realistic atomic models to describe these ef- 
fects, especially for atoms in a specific external poten- 
tial giving rise to a Coulomb instability. To this end we 
have propos ed a generalized non-rigid shell model (see 
Refs l7ll2l2^ . The model represents a variational method 
for the many-electron atomic configurations with the trial 
parameters being the coordinates of the center of the one- 
particle atomic orbital. The resulting displacement of 
the atomic orbitals allows a simple interpretation of the 
electron density redistribution stemmed from taking into 
account the electron-electron repulsion, and the symme- 
try of a system can be readily used for the construction 
of the trial many-electron wave function. As a whole, 
the model bears a strong resemblance to the conventional 
well-known shell model by Dick and Overhauser 26 widely 
used in lattice dynamics. In frames of the model the ionic 
configuration with filled electron shells is considered to 
be constituted of an outer spherical shell of 2(21+1) elec- 
trons and a core consisting of the nucleus and the remain- 
ing electrons. In an electric field the rigid shell retains 
its spherical charge distribution but moves bodily with 
respect to the core. The polarizability is made finite by 
a harmonic restoring force of spring constant k which 
acts between the core and shell. The shells of two ions 
repel one another and tend to become displaced with re- 
spect to the ion cores because of this repulsion. The 
respective displacement vector appears to be a simplest 
collective coordinate which specifies the change of the 
electron-nucleus attraction. It should be noted that such 
a displacement does not imply any variation in electron- 
electron repulsion and respective correlation energy. 

However, as we shall see below, a simple shell model 
can be easily generalized to take account of correlation 
effects. To this end we must consider the displacements 
of separate one-electron orbitals to form the set of the 
variational parameters in a correlation function. Then 
we can introduce both the displacement of the center 
of " gravity" for filled shell and a set of the relative dis- 
placements of separate one-electron orbitals with regard 
to each other. The former form an "acoustical" mode and 
are described in frames of conventional shell model, while 
the latter form different novel "optical" modes. Such a 
seemingly naive non-rigid shell picture can provide both 
the microscopic substantiation of the conventional shell 
model and its generalization. Moreover, this non-rigid 
shell model points to a physically clear procedure to ac- 
count for the correlation effects. Indeed, the "optical" 
displacement mode is believed to provide the minimal 



electron-electron repulsion. The non-rigid shell repre- 
sents a novel specific atomic state that can be remarkably 
realized near a Coulomb instability. 

The idea of displaced atomic shells has appeared many 
years ago^i in the very beginning of the quantum chem- 
ical era, and reflected the naive picture of the repelling 
electrons. However, the physically sound idea did not re- 
ceive the relevant position in the hierarchy of correlation 
effects. 



2. Electron-lattice relaxation effects 

As it is mentioned above, the minimal energy cost of 
the optically excited disproportionation or electron-hole 
formation in insulating cuprates is 2.0 — 2.5 eV. However, 
the question arises, what is the energy cost for the ther- 
mal excitation of such a local disproportionation? The 
answer implies first of all the knowledge of relaxation en- 
ergy, or the energy gain due to the lattice polarization 
by the localized charges. The full polarization energy 
R includes the cumulative effect of electronic and ionic 
terms, associated with the displacement of electron shells 
and ionic cores, respectively^ The former term R op t is 
due to the non-retarded effect of the electronic polar- 
ization by the momentarily localized electron-hole pair 
given the ionic cores fixed at their perfect crystal posi- 
tions. Such a situation is typical for lattice response ac- 
companying the Franck-Condon transitions (optical ex- 
citation, photoionization) . On the other hand, all the 
long-lived excitations, i.e., all the intrinsic thermally ac- 
tivated states and the extrinsic particles produced as a 
result of doping, injection or optical pumping should be 
regarded as stationary states of a system with a deformed 
lattice structure. These relaxed states should be deter- 
mined from the condition that the system energy has a 
local minimum when account is taken of the interaction 
of the electrons and holes with the lattice deformations. 
At least, it means that we cannot, strictly speaking, make 
use of the same energy parameters to describe the optical 
(e.g. photoexcited) hole and thermal (e.g. doped) hole. 

For the illustration of polarization effects in cuprates 
we apply the shell model calculations to look specifically 
at energies associated with the localized holes of Cu 3+ 
and in "parent" La2Cu04 compound. It follows from 
these calculations that there is a large difference in the 
lattice relaxation energies for O - and Cu 3+ holes. The 
lattice relaxation energy, -AR" h , caused by the hole local- 
ization at the O-site (4.44 eV) appears to be significantly 
larger than that for the hole localized at the Cu-site (2.20 
cV) . This indicates the strong electron-lattice interaction 
in the case of the hole localized at the O-site and could 
suggest that the hole trapping is more preferential in the 
oxygen sublattice. In both cases we deal with the several 
cV-cffcct both for electronic and ionic contributions to 
relaxation energy. Moreover, such an estimation seems 
to be typical for different insulators. 4,5 It is worth not- 
ing that the electron-lattice interaction is believed to be 
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one of the main sources of correlated particle hopping re- 
sulting in different probabilities for two types of a charge 
transfer. 



3. Generalized Peierls- Hubbard model and "negative-U" 
effect 

Transition metal oxides with strong electron and lat- 
tice polarization effects need in a revisit of many conven- 
tional theoretical concepts and approaches. In particu- 
lar, we should modify conventional Hubbard model as 
it is done, for instance, in a "dynamic" Hubbard model 
by Hirsch— or a modified Peierls-Hubbard models with 
a classical description of the anharmonic core/shell dis- 
placements. Having in mind the application to insulating 
cuprates let address a square lattice Hubbard model with 
a half-filling and a strong on-site coupling of valent hole 
with core/shell displacements, which is described by the 
following Hamiltonian 

H = - ]T t{ij)c\ a c ja + U^2 n i° n iv' 

<i-j>,c i.aa' 



+ '^ J Van(qi,ni) + ^ v int(qi,qj), (8) 
i <i,j> 

where c\ a (cj a ) are creation (annihilation) operators for 
valent hole; t(ij) = t(qi,qj) is the transfer integral be- 
tween two neighboring lattice sites which depends on the 
dimensionless core/shell displacement coordinate; U is 
the on-site repulsion energy; v an (qi,rii) is the configu- 
rational energy that incorporates the coupling between 
valent holes and the site-localized anharmonic core/shell 
mode with dimensionless displacement coordinate g^; 

Van(qi,ni) = a{ni)ql - b(rii)qf + c(n i )g l 6 , (9) 

where a, b, c are the functions of the hole occupation num- 
ber such as 

a(rii) = a + ain. t + a 2 nf , (10) 

It is clear that v an (qi,rii) includes the renormalization 
both of the one-particle energy and the on-site hole-hole 
repulsion. The last term in i)22|) represents the intersite 
configurational coupling. The q-dependence of transfer 
integral implies the correlated character of the hole hop- 
ping, and can be transformed into the effective depen- 
dence on hole occupation number— 

t(rii, rij) = t(l + a(rii + rij) + flniTij). (11) 

with a, (3 being the correlated hopping parameters. 

The conventional Hubbard Hamiltonian, or t-U-mode\, 
stabilizes the spin density wave (SDW) electron order 
with m = 1. In a strongly correlated limit U > t the 



Hubbard model reduces to a Heisenberg antiferromag- 
netic model. Depending on the parameters of the hole- 
configurational coupling and correlated hopping the mod- 
ified Hubbard Hamiltonian (1221) can stabilize the "dis- 
proportionated" or charge ordered (CO) electron phase 
with the on-site filling numbers n = 0, and n = 2 thus 
leading to the "negative-U" effect. Even simple modified 
model turns out to be very complicated and leads to a 
very rich physics. 23 Depending on the values of parame- 
ters the system yields the SDW phase with no core/shell 
displacements as a true ground state with a global min- 
imum of free energy, and CO phase with shell displace- 
ments as a false ground state with a local minimum, or 
vice versa£ Strong anharmonicity v an (qi, m) makes pos- 
sible phase transitions between the phases the first order 
ones. 



4- Vibronic reduction of charge transfer integrals 

In general, charge, spin and vibronic modes are 
strongly coupled and so we have to do with the hybrid 
modes. For a weak intermode coupling regime the charge 
transfer is accompanied by the induced local structural 
fluctuations, that provides the vibronic reduction of the 
charge transfer integral: 

ti2=*S ) K vib K mb = ( X i|X2) 2 , (12) 

where K v n is a vibronic reduction factor, (X1IX2) is an 
overlap integral for the local oscillatory states with and 
without particle transferred. In an opposite regime of 
the strong intermode coupling one assumes that different 
electronic parameters for the e- and /i-centers are distin- 
guished significantly up to different type of the adiabatic 
potential and appropriate JT mode. This regime favors 
the charge localization. 

The vibronic reduction factor K v a,, or the Franck- 
Condon factor— may be written as follows 

K mb =N exp(--y), (13) 

where N and 7 in a complicated manner depend on 
the vibronic constants, the oxygen and 3d-metal atomic 
masses. For a simplest one-dimensional single-mode case 

*— j^-pf-j^). (») 

where l\ and I2 are the effective oscillatory lengths of 
the 1- and 2-centers, respectively, r = h/1%, AQ is the 
distance separating the minima of the adiabatic potential 
for the 1- and 2-centers. 



5. Spin reduction of charge transfer integrals 

Overall the paper we neglect a spin degree of freedom 
which can crucially impact on the charge transport. The 
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most part of 3d oxides are characterized by an antiferro- 
magnetic spin background that implies a localization ef- 
fect due to strong spin reduction of one-particle transfer 
integrals and the probability amplitude for a polar cen- 
ter transfer M ± + M° -> M° + M ± or the motion of the 
electron (hole) center in the matrix of M°-centers. From 
the other hand, in antiferromagnets there is no problems 
with another type of the CT which specifies the probabil- 
ity amplitude for a local disproportionation, or the spin- 
singlet e/i-pair creation: M° + M° -> M ± + M T , and 
the inverse process of the spin-singlet e/i-pair recombina- 
tion: M ± + M T -> M° + M°. In other words, the spin 
subsystem can strongly affect the correlated character of 
the charge transfer leading to unconventional situations 
like that of spin-singlet eh-pairs moving through the lat- 
tice freely without disturbing the antiferromagnetic spin 
background, in contrast to the single particle motion. So, 
it seems that the situation in antiferromagnetic 3d insu- 
lators may differ substantially from that in usual semi- 
conductors or in other bandlike insulators where, as a 
rule, the effective mass of the electron-hole pair is larger 
than that of an unbound electron and hole. 



III. S=l PSEUDOSPIN FORMALISM FOR 
MODEL MIXED VALENCE SYSTEM 

A. Pseudospin operators 

The problem of the multi-stability of solids looks rather 
trivial when one say about the orbital and/or spin degrees 
of freedom. Usually in such a case we start from the lat- 
tice of coupled orbital and/or spin momenta described 

I 



Here we make use of standard symbols for the Wigner 
coefficients and reduced matrix elements. In a more 
conventional Cartesian scheme a complete set of the 
non-trivial pseudo-spin operators would include both S 
and a number of symmetrized bilinear forms {SiSj} = 
(SiSj + SjSi), or spin-quadrupole operators, which are 
linearly coupled to and V 2 , respectively 

Vq = S q ] So = S z , S± = ^-j=(S x ± iS y ) : 

V 2 ex (3Sf-S 2 ),V& cx (S z S ± +S ± S z ),Vl 2 ex S 2 ± . (16) 

These pseudo-spin operators are not to be confused with 
real physical spin-operators; they act in an imaginary 
pseudo-space. 



by the relevant (spin-)Hamiltonian that implies the va- 
riety of possible collective orbital and/or spin orderings 
that compete with each other under different external 
conditions. In other words, the multi-stability accom- 
panies the basic degeneracy inherent to a certain atom, 
ion, or center with a nonzero orbital and/or spin momen- 
tum. Such an outlook is believed to be easily extended 
to systems with charge degree of freedom which can be 
represented to be a system of either centers which pos- 
sible charge states form a pseudo-multiplet. Below we 
address a simple model of a mixed-valence system with 
three possible stable valent states of a cation-anionic clus- 
ter, hereafter M: M°,M ± , forming the charge (isospin) 
triplet. Starting from M° state as a bare vacuum state, 
we may address the M ± centers as a result of pseudo- 
spin AS Z = ±1 deviation, or as a hole and electron, 
respectively. Below we intend to concentrate themselves 
on charge degree of freedom, and that is why we neglect 
the orbital, spin, and lattice degrees of freedom. It im- 
plies a renormalization of different parameters, mainly it 
concerns the charge transfer. 

Similarly to the neutral-to-ionic electronic-structural 
transformation in organic charge-transfer crystals (see 
paper by T. Luty in RefllH) the system of charge triplets 
can be described in frames of the S=l pseudo-spin for- 
malism. To this end we associate three charge states of 
the M-center with different valence: M°,M with three 
components of S = 1 pseudo-spin (isospin) triplet with 
Ms = 0, +1, —1, respectively. 

The S = 1 spin algebra includes three independent 
irreducible tensors of rank k = 0, 1, 2 with one, three, 
and five components, respectively, obeying the Wigner- 
Eckart theorem^ 



-Ms k q M' s )^ k WS). d5) 
I 

To describe different types of pseudo-spin ordering in 
a mixed-valence system we have to introduce eight or- 
der parameters: two diagonal order parameters (S z ) 
and (S z ), and six off-diagonal order parameters (V^) 
(q =/= 0). Two former order parameters can be termed 
as valence and ionicity, respectively. The off-diagonal 
order parameters describe different types of the valence 
mixing. Indeed, operators Vg (q ^ 0) change the z- 
projection of pseudo-spin and transform the \SM$) state 
into \SMs + q) one. In other words, these can change 
valence and ionicity. It should be noted that for the 
5=1 pseudospin algebra there are two operators: V±± 
and V±i, that change the pseudo-spin projection by ±1, 
with slightly different properties 

(0|£±|=Fl> = (±l|&fc|0)==Fl, (17) 



(SM s \V q k \SM s ) 



-1) 



S-M S 
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but 



(0\(S z S ± +S ± S z )\Tl) = -(±l\(S z S ± + S ± S z )\0) =+1. 

(18) 



B. Gell-Mann operators and generalized 
pseudospin Hamiltonian 

Three spin-linear (dipole) operators Si, 2,3 and five in- 
dependent spin-quadrupole operators {Si,Sj} — fS 2 ^ 
given S = 1 form eight Gell-Mann operators being the 
generators of the SU(3) group. Below we will make use of 
the appropriate Gell-Mann 3x3 matrices A( fc ), which dif- 
fer from the conventional only by a remuneration; 30 
A(D = a< 6 >, A« = A( 3 ), A< 3 > '= AW, A< 4 ) = A( 5 \ 
A (5) = _ A (2), AM = AW, A( y ) = A«, AW = A«. 
First three matrices A^ 1 ' 2,3 ) correspond to linear (dipole) 
spin operators: 



A^=S y : A^=S Z 



while other five matrices correspond to quadratic 
(quadrupole) spin operators: 

A^ = -{S z S y }; A^ = ~{S X S Z }; A< 6 > = -{S x S y }; 



A (7) = -^{Sl 



si - 2Sf 



A (8) 



2£ 



with B being a unit 3x3 matrix. 

The generalized spin-1 model can be described by the 
Hamiltonian bilinear on the SU(3)-generators A^ 



H = - 



i,7] k,m=l 



(k) 



(19) 



Here i,r/ denote lattice sites and nearest neighbors, re- 
spectively. This is a S = 1 counterpart of the 5=1/2 
model Heisenberg Hamiltonian with three generators of 
the SU(2) group or Pauli matrices included instead of 
eight Gell-Mann matrices. 



±|1 — 1 >)/V2. The linear (dipole) pseudo-spin operator 
within \x, y, z) basis is represented by a simple matrix: 

< 4>i\Sj\ip k >= -iEij k , 

and for the order parameters one easily obtains: 

< S >= -2[a x b]; < >= 2(% - a^j - hbj) 

(21) 

given the normalization constraint a 2 -I- b 2 = 1. Thus, 
for the case of spin-1 system the order parameters are de- 
termined by two classical vectors (two real components 
of one complex vector c = a + ib from l|2U|l ) . The two 
vectors are coupled, so the minimal number of dynamic 
variables describing the S = 1 spin system appears to 
be equal to four. Along with a, b vectors one might 
introduce 1 = [a x b] . Hereafter we would like to empha- 
size the director nature of the c vector field: ip(c) and 
c) describe the physically identical states. It is worth 
noting that the coherent states provide the optimal way 
both to a correct mean-field approximation (MFA) and 
respective continuous models^ 
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FIG. 2: Schematic energy spectrum of two nn M-centers sys- 
tem (see text for details). Arrows mark the dipole-allowed 
CT transitions. 



IV. MODEL MIXED- VALENCE SYSTEM 



C. Generalized mean- field model 

In frames of a classical, or mean-field description of 
the S = 1 quantum pseudo-spin system we start from a 
coherent state approximation with trial functions^. 

V> = II *C?')ifc = II (20) 

j£lattice jdattice 

Here j labels a lattice site and the spin functions ipi in a 
Cartesian basis are used: ip z = j 10 > and i/) Xy y ~ (|11 > 



A. Effective pseudo-spin Hamiltonian 

Effective pseudo-spin Hamiltonian for our model 
mixed-valence system should incorporate a large num- 
ber of contributions that describe different long- and 
short-range coupling between M ' centers, single-ion 
and two- ion terms. Single-site terms can be subdivided 
into single-ion anisotropy and pseudo-Zeeman interac- 
tion. Bilinear and biquadratic two-site terms can be sub- 
divided into diagonal interactions like " density-density" , 
and off-diagonal terms that describe charge fluctuations 
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conserving the total charge of the system, such as one- 
electron(hole) and two-electron(hole) transport. An ef- 
fective pseudo-spin Hamiltonian of the model mixed- 
valence system which takes into account the main part of 
aforementioned contributions can be represented as fol- 
lows 

H^iAiSl-hiS^+Y, VijSlSl+Y, V l3 S tz S JZ + 

i i .j ■ ■ i .J ■ 

[D§\s i+ Sj- + SiS j+ ) + Dff(T i+ Tj_ + T t _T J+ )] 

<i,j> 

+ tij(S? + S]_ + Sf_S] + ), (22) 
<i,3> 

where 

T± = (S Z S± + S±Sg). 

Two first single-ion terms describe the effects of bare 
pseudo-spin splitting, or the local energy of M 0,± cen- 
ters. Interestingly, the parameter A can be related with 
correlation Hubbard parameter U: U — 2A. The sec- 
ond term may be associated with an external, generally 
speaking, non- uniform pseudo-magnetic field hi, in par- 
ticular, a real electric field. It is easy to see that it de- 
scribes an electron/hole assymetry. The third and fourth 
terms describe the effects of long- and short-range inter- 
ionic interaction including screened Coulomb and cova- 
lent coupling. 

If to apply the familiar spin terminology, the first 
term in (|22|l represents a single-ion anisotropy, the sec- 
ond does the Zeeman term, the fourth and fifth do the 
anisotropic Heisenberg exchange, and the third and sixth 
do the biquadratic spin-quadrupolar coupling. To il- 
lustrate the role of different terms in (|22|l we present 
in Fig.l a schematic energy spectrum of nn pair of 
M centers provided an eh-symmetry (h — 0) and 1 00) 
ground state (A > 0). It is worth noting the effect of a 
renormalization of the ground state due to eh-pair cre- 
ation/recombination effect (ty 0) with a stabilization 
energy S « \t'^\ 2 /2A. Two electron-hole states with S- 
(even) and P- (odd) type symmetry have a very strong 
dipole coupling with the large value of S — P transition 
dipole matrix element: 

d=|(S , |d|P)|«2e J R M M«2ex4A. (23) 

Contrary to P-type pair state the S-type one is dipole- 
forbidden and corresponds to a so-called two-photon 
state. 

One should note that despite many simplifications, and 
first the neglect of orbital and spin degrees of freedom, 
quenched lattice approximation, the effective Hamilto- 
nian (|22|l is rather complex, and represents one of the 
most general forms of the anisotropic S = 1 non- 
Heisenberg Hamiltonians. For the system there are two 



classical (diagonal) order parameters: (S z ) = n being 
a valence, or charge density with electro-neutrality con- 
straint J2i n i — J2i = 0' an d (S z ) = n p being the 
density of polar centers M , or "ionicity". In addition, 
there are two unconventional off-diagonal order param- 
eters: "fermionic" (5+) and "bosonic "(S+); the for- 
mer describes a phase ordering for the disproportiona- 
tion reaction, or the single-particle transfer, while the 
latter does for exchange reaction, or for the two-particle 
transfer. Indeed, the S + operator creates a hole and is 
fermionic in nature, whereas the does a hole pair, and 
is bosonic in nature. 



B. Single and two-particle transport 

The last three terms in (I22fl representing the one- and 
two-particle hopping, respectively, are of primary impor- 
tance for the transport properties, and deserve special 
interest. 

Two types of one-particle hopping are governed by two 
transfer integrals D^' 2 \ respectively. In accordance with 
(|T7|) and (PJ) the transfer integral t'^ = (D^f + D^f) 
specifies the probability amplitude for a local dispropor- 
tionation, or the eh-pair creation 

M°+M° -> M ± +M T ; 

and the inverse process of the eh-pair recombination 

M ± + M T -> M° + M°, 

while the transfer integral t^ = (D^ — D^) specifies 
the probability amplitude for a polar center transfer 

M ± + M° — » M° + M ± , 

or the motion of the electron (hole ) center in the matrix 
of M° -centers or motion of the M°-center in the matrix 
of M ± -centers. It should be noted that, if t'^ — but 
tfj 7^ 0, the eh-pair is locked in two-site configuration. 

The two-electron(hole) hopping is governed by transfer 
integral tij, or a probability amplitude for the exchange 
reaction: 

M ± +M T -> M =F + M ± . 

or the motion of the electron (hole ) center in the matrix 
of hole (electron) centers. 

It is worth noting that in Hubbard-likc models all the 
types of one-electron(hole) transport are governed by 
the same transfer integral: t'^ = if!- = t^, while our 
model implies independent parameters for a dispropor- 
tionation/recombination process and simple quasiparti- 
cle motion in the matrix of Af°-centers. In other words, 
we deal with a "correlated" single particle transport £ 
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C. Mean-field approximation:three generic 
MFA-phases 

First of all we would like to emphasize the differ- 
ence between classical and quantum mixed- valence sys- 
tems. Classical (or chemical) description implies the 
neglect of the off-diagonal purely quantum CT effects: 
Di 1 ' 2 ) — t = 0, hence the valence of any site remains to 
be definite: 0, ±1, and we deal with a system of localized 
polar centers. In quantum systems with a nonzero charge 
transfer we arrive at quantum superpositions of different 
valence states resulting in indefinite on-site valence and 
ionicity which effective, or mean values (S z ) and (S 2 .) can 
vary from —1 to +1 and to +1, respectively. 

Making projection of the effective pseudo-spin Hamil- 
tonian for the system onto a space of states like l|20() , we 
obtain an energy functional which equivalent to a clas- 
sical energy of the two coupled vector (a, b) fields de- 
fined on the common lattice. Thus, in the framework of 
the pseudo-spin 5 = 1 centers model when the collective 
wave function is represented to be a product of the site 
functions like (J2DJ, the quantum problem is reduced to a 
classical variation problem for a minimum of the energy 
for two coupled vector fields. 

In frames of mean-field approximation (MFA) we may 
make use of coherent states (PU)l that provide a physically 
clear assignment of different phases with a straightfor- 
ward recipe of its qualitative and quantitative analysis. 
All the MFA-phases one may subdivide into those with 
a definite and indefinite ionicity, respectively. There are 
two MFA-phases with definite ionicity; 

1) Insulating monovalent M°-phase with (S 2 ) = 

0: 

The M°-phase is specified by a simple uniform arrange- 
ment of a and b vectors parallel to z-axis: a || b || O z . 
In such a case the on-site wave function is specified by 
unit vector (a, or b) parallel to z-axis. It is a rather 
conventional ground state phase for various charge trans- 
fer insulators such as oxides with a positive magnitude 
of A parameter (U > 0). All the centers have the 
same bare M° valence state. In other words, the M°- 
phase is characterized both by definite site ionicity and 
valence. So, all the order parameters turn into zero: 
(S z ) = (S 2 ) = (S+) = (Sl) = 0. This is an "easy-plane" 
phase for the pseudo-spins, but an "easy-axis" one for 
the a and/or b vectors. 

2) Mixed-valence binary (disproportionated) 
M ± -phase with (S 2 ) = 1: 

This phase usually implies an overall disproportiona- 
tion M° + M° -> M ± + M T that seems to be realizable 
if A parameter becomes negative one (negative U < 
effect). It is a rather unconventional phase for insula- 
tors. All the centers have the "ionized" valence state, 
one half the M + state, and another half the M~ one, 
though one may in common conceive of deviation from 
fifty-fifty distribution. A simplified "chemical" approach 
to Af ± -phase as to a classical disproportionated phase 
is widely spread in solid state chemistry. 14 In contrast 



with the M° phase the M ± -phase is specified by a pla- 
nar orientation of a and b vectors (a, b 1 2 ) with a 
varied angle in between. There is no fermionic trans- 
port: (5+) = 0, while the bosonic one may exist, and, 
in common, (S£) = - cos(0 a - <fo)e^o + ^ ^ 0. This is 
an "easy-axis" phase for the pseudo-spins, but an "easy- 
plane" one for the a and b vectors. 

The mixed valence phase as a system of strongly 
correlated electron and hole centers is equivalent to the 
lattice hard-core Bose system with an inter-site repulsion, 
or electron-hole Bose liquid (EHBL) in contrast with EH 
liquid in conventional semiconductors like Ge, Si where 
we deal with a two-component Fermi- liquid. Indeed, one 
may address the electron M~ center to be a system of a 
local boson (e 2 ) localized on the hole M + center: M ~ = 
M+ + e 2 . 

In accordance with this analogy we assign three well 
known molecular-field uniform phase states of the 
binary mixture: 

i) charge ordered (CO) insulating state with 
(S z ) = ±1, a _L b, and zero modulus of bosonic off- 
diagonal order parameter: | (>S'^_ ) | = 0; 

ii) Bose-superfluid (BS) superconducting state 
with (S z ) =0, a and b being collinear, = e 21 *; 

iii) mixed Bose-superfluid-charge ordering 
(BS+CO) superconducting state (supersolid) 
with < KS^)! < 1, a and b being oriented in xy-plane, 
but not collinear, (S%) = - cos(0 a - ^e^ 4 ^ ^ 0. 

In addition, we should mention the high-temperature 
non-ordered (NO) Bose-metallic phase with <C S z ^>= 0. 

Rich phase diagram of binary mixture with uncon- 
ventional superfluid and supersolid regions looks tempt- 
ing, however, actually, their stabilization requires strong 
suppression of Coulomb repulsion between electron (hole) 
centers. Despite significant screening effect, the stabi- 
lization of uniform BS or BS+CO superconducting state 
as a result of a disproportionation reaction in a bare 
insulator— seems to be unrealistic. 

3) Mixed- valence ternary ("under- 
disproportionated" ) M 0,± -phase: 

For two preceding cases the order parameter (S 2 ), 
or ionicity relates to its limiting values (0 or 1, re- 
spectively). For the MFA-phase with indefinite ionicity, 
or mixed- valence ternary ( "under-disproportionated" ) 
M 0,± -phase, < (S 2 ) < 1, that is we have a mixture 
of the M°, centers. From the viewpoint of the clas- 
sical a, b vectors formalism the phase corresponds to the 
arbitrarily space-oriented 1 = [a x b] vector. Both off- 
diagonal order parameters, fermionic (S+) and bosonic 
(S 2 .) are, in common, non-zero, albeit with some correla- 
tion in between. So, for the ternary system one expects 
a coexistence of fermionic and bosonic transport. 

It should be noted that a complete pseudo-spin descrip- 
tion of the two last model mixed- valence systems implies 
a two-sublattice approximation to be a minimal model 
compatible with a sign of Coulomb interaction and a re- 
spective tendency towards the checkerboard-like charge 
ordering. 
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V. INSULATING MONOVALENT M°-PHASE 

Insulating monovalent M°-phase is a typical one for 
the ground state of insulating transition metal oxides, 
or Mott- Hubbard insulators. It is worth noting that 
in frame of conventional band model approach the M°- 
phase, e.g., in parent cuprates, is associated with a metal- 
lic half-filled hole band system. Below we address differ- 
ent types of quasiparticle excitations in such a system fo- 
cusing on the features of the correlated hopping, governed 
generally by two transfer integrals and 1/L which com- 
petition results in unconventional properties of electrons 
and holes in bare insulating monovalent M°-phase. We 
show that M°-phase can be unstable with regard to the 
self-trapping of CT excitons and nucleation of droplets 
of EH Bose liquid. 

A. Electron, hole, and electron-hole excitations 

Starting from monovalent M°-phase as a vacuum state 
|0) we introduce an electron- hole representation where 
Mr center is derived as a result of an electron creation 

df |0), and center is derived as a result of a hole cre- 
ation b\ |0) . Then we transform pseudo-spin Hamiltonian 
(|22|l into that of a system of effective electrons and holes 

H = Aj2(n1 + <) + E *« ( fi fo + 

E W** W ) ( n i n i + (ij'K»i + V Ch (ij) (»i »i + n)nX)\ 
hi 

+ E4( a !^ + a A-) (24) 
(ij) 

where Vu — > oo to prohibit two-particle occupation of 
a single site. Here we suppose h = that provides an 
electron-hole symmetry. The first line in 124|) represents 
the single particle (electrons/holes) terms, the second one 
does the interparticle coupling, while the third one de- 
scribes the creation and annihilation (recombination) of 
ch-pairs. It is worth noting that the latter terms describe 
some sort of eh-coupling. 

In terms of a pseudo-spin analogy the electrons and 
holes are associated with pseudo-spin AS Z — ±1 devia- 
tions for an easy-plane magnet, localized or delocalized 
(pseudo-spin wave). 

The behavior of electron/hole system crucially depends 
on the relation between two transfer integrals t'^ , t"j Be- 
low we address two distinct limiting situations: 

I. t"j = 0: Forbidden recombination/creation regime. 
In this case we deal with the bands of well defined elec- 
trons and holes with a charge gap Eg> h = A — z\t' nn \ for 
both types of carriers. Optical gap for unbound eh-pairs 
is Eg h = 2Eg' h . However, in such a case we may expect 



for the formation of Wannier excitons, or eh-pairs bound 
due to a screened Coulomb eh-coupling. In terms of 
pseudo-spin formalism the Wannier excitons may be re- 
garded as two pseudo-spin waves having formed a quasilo- 
calized state due to a long-range antiferromagnctic Ising 
exchange: V i3 S iz S jz . 

II. — 0: Regime of localized electrons and holes 
with a dimerization effect and well defined nn eh-pairs, 
or CT excitons. In such a case the charge gap is Eg- h = A 
for both types of localized quasiparticles. As it is clearly 
seen in Fig.l the t"j ^ hopping results in a dimeriza- 
tion effect with a quantum renormalization of the vacuum 
state and indefinite ionicity, the formation of two types 
of localized eh-pairs, or CT excitons of S- and P-type. 
In frames of our nn approximation the CT excitons are 
localized. Optical gap is determined by the energy of 
P-type CT exciton: Ef = 2A + S. 

Thus we arrive at two limiting types of monovalent 
M insulators with a dramatic difference in behavior of 
electrons and holes, as well as electron-hole pairs. In type 
I insulators (tL ^> t"j) we deal with well defined bands 
of electrons and holes forming Wannier excitons, while 
in type II insulators (t"j ^> t' io ) we deal with localized 
electrons and holes which can form nn eh-pairs, or CT 
excitons. 

The most part of 3d oxides are characterized by an 
antiferromagnetic spin background that implies strong 
spin reduction of one-particle transfer integrals t'^. In 
other words, these, seemingly, belong to type II insula- 
tors, where spin-singlet CT excitons can move through 
the lattice freely without disturbing the antiferromag- 
netic spin background, in contrast to the single hole mo- 
tion. So, it seems that the situation in antiferromagnetic 
3d insulators differs substantially from that in usual semi- 
conductors or in other bandlike insulators where, as a 
rule, the effective mass of the electron-hole pair is larger 
than that of an unbound electron and hole. 

The Wannier excitons are formed due to an ch- 
coupling, while the CT excitons are formed due to a ki- 
netic cutoff, or a specific feature of correlated hopping, 
in other words, the former have a potential, while the 
latter a kinetic nature. 

It is worth noting that both M-centers within P-type 
CT excitons have a certain ionicity in contrast to S-type 
CT excitons which can mix with bare M°M° ground 
state. CT excitons form peculiar quanta of a dispropor- 
tionation reaction and may be viewed to be a minimal 
droplet of electron-hole Bose liquid. 

In general, eh-excitations in M° phase consist of su- 
perpositions of pairs of free electrons and holes, and CT 
excitons. One expects two types of superpositions: CT 
exciton-like and band-like. The former have a localized 
character, while the latter have an itinerant one. 

The nature of the optical excitations accompanied by 
creation of electron-hole pairs in 3d oxides is not fully 
understood. One of the central issues in the analysis 
of electron-hole excitations is whether low-lying states 
are comprised of free charge carriers or excitons. A con- 
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ventional approach implies that if the Coulomb interac- 
tion is effectively screened and weak, then the electrons 
and holes are only weakly bound and move essentially 
independently as free charge-carriers. However, if the 
Coulomb interaction between electron and hole is strong, 
excitons are believed to form, i.e. bound particle-hole 
pairs with strong correlation of their mutual motion. 

One of the most popular criteria to discriminate be- 
tween the states relates to the band gap: states below 
the charge gap correspond to excitons with binding en- 
ergy Eb — E g ~ E, and states above the charge gap do to 
free electron-hole pairs. However, this criteria seems to 
be oversimplified, and the states should be characterized 
as bound or unbound according to the scaling of the av- 
erage electron-hole separation with system size. Excitons 
are entities with small electron-hole separation which re- 
main finite as the system size is increased. By contrast, 
the average separation between two free charge carriers 
increases indefinitely with system size. To distinguish 
bound and unbound electron-hole states one might use 
the density-density correlation function^ 

C(i,j) = ((hi - (hi^hj - (hj))}, 

which measures a correlation of charge fluctuations on 
site i to a charge fluctuations on site j. A negative value 
correlates an excess (deficit) of charge with deficit (ex- 
cess), or electron- hole distribution. In frame of pseudo- 
spin approach this correlation function measures the lon- 
gitudinal (|| z) short-range antiferromagnetic fluctuations 

C(i,j) = (S iz S jz )-(S iz )(S jz ). 

The results of the direct computations for conjugated 
polymers in frames of the extended Hubbard model^i 
show for different symmetry sectors that there exist un- 
bound states below the charge gap, and bound states 
above the charge gap. Thus, the charge gap, often used 
to define the binding energy of excitons, is not a decisive 
criterion by which to decide whether a state is bound or 
unbound. 

In practice, many authors consider excitons to con- 
sist of real-space configurations with electrons and holes 
occupying the nearest neighbor sites, while the elec- 
trons and holes are separated from each other in the 
conduction-band states^. 

The relative energy position and transition intensity of 
excitons and free unbound electron-hole pairs is an issue 
of large complexity. For instance, in polydiacetylenes^ 
the large absorption peak at 1.85 eV is attributed to an 
exciton, since photoconductivity is absent in this energy 
region. The latter has a threshold at around 2.4 eV. In- 
terestingly, that observed photoconductivity band is not 
visible in the conventional linear (one-photon) absorp- 
tion spectrum. Because of the strong oscillator strength 
of the exciton, the conduction band has a weak oscillator 
strength and is enveloped by the high-energy tail of the 
exciton peak. 



B. Coupling with electromagnetic field 

The electron-hole excitations and optical properties of 
strongly correlated electron systems is of current both 
experimental and theoretical interest. In particular, the 
optical conductivity is of fundamental theoretical interest 
because the spectral weight at low frequencies seems to 
be the natural order parameter for the Mott transition. 
The great body of experimental information regarding 
the eh-excitations in solids is provided by optical and 
electron energy loss (EELS) spectroscopy, however, its in- 
terpretation depends crucially on the theoretical scheme 
used. The optical conductivity of different model systems 
has been studied by approximate mean field calculations, 
by analysis of integrable ID models, by exact diagonal- 
ization of small systems, and by quantum Monte Carlo 
techniques. The uncertain quantitative applicability of 
the analytic mean field calculations, the size limitations 
of the exact diagonalization and Monte Carlo results and 
the possibility that the integrable ID models do not ex- 
hibit generic behavior one lead to consider other methods 
for obtaining information about the optical conductivity. 
An outgoing beyond the effective Hamiltonian methods 
with restricted basis which are relevant for description 
of the lowest energy excitations in the extremely limited 
range of energies, and the elaboration of effective meth- 
ods to describe optically excited states is a challenging 
problem in solid state physics. 

Making use of a standard Peierls transformation in 
hopping terms in (|24|) we arrive at an effective Hamil- 
tonian for the coupling with electromagnetic field 

i.ij ► t^jc J , (25) 

- = -± £ 3 A(r)dr, (26) 

where A is the vector-potential, and integration runs over 
line binding the i and j sites. For the nearest neighbours 
one may use the simplified relation 

($ J --$ i ) = -^(A(R i )"Rji)- (27) 

The current density operator one may represent to be 
a sum of three terms 

j ee (R,) = ^| *ii*MOifij - fiifi}). (28) 
3 

j hh (Ki) ttfRtf(&& - kb% (29) 

3 



fCRO = ^ £ UJUj{H\b\ - k&i), (30) 

3 
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where the two first terms describe the electron and hole 
currents, respectively, while the third term describe the 
current fluctuations due to eh-pair creation/annihilation. 

Standard linear-response theory then yields for optical 
conductivity (T = 0) 

Rea(q,^) = (tt/w) £ \^ ge {a L )\ 2 5{E e - E g - foj), (31) 

e 

where 

M S e(q) = (*gbq|*e) 

is a transition matrix element, and sum runs over all the 
excited ^ e states. 

The first two current density operators Q28 [l -I|29 [l de- 
scribe electron and hole intraband transition so that the 
optical absorption spectra in M° phase are specified only 
by the latter term (|30|l . In other words, the optical re- 
sponse in M° phase of our model system is determined 
only through a CT exciton channel. 

It should be noted that the electron/hole current op- 
erators in can t> e expressed through a superpo- 
sition of antisymmetric pseudospin operators: 

(Si+Sj- - Si.S j+ ) - TUT,- - f i _f j+ ), (32) 

while eh-current fluctuation operator in (|30|l can be ex- 
pressed through another superposition of the same anti- 
symmetric pseudospin operators 

(Si+Sj- — Si~Sj+) + Ti + Tj- — Ti-Tj + ). (33) 

C. Charge transfer instability and CT exciton 
self-trapping 

Electron and hole in a CT exciton in type II M° insu- 
lator are strongly coupled both in between and with the 
lattice. In contrast with conventional wide-band semi- 
conductors where the excitons dissociate easily produc- 
ing two-component electron-hole gas or plasma^ small 
CT excitons both free and self-trapped are likely to be 
stable with regard the eh-dissociation. To illustrate the 
principal features of CT exciton self-trapping effect we 
address a simplified two-level model of a two-center MM 
cluster in which a ground state and a CT state are as- 
sociated with a pseudospin 1/2 doublet, | f) and | |), 
respectively. In addition, we introduce some configura- 
tional coordinate Q associated with a deformation of the 
cluster, or respective anionic background. 7 Such a model 
is typical one for so-called (pseudo)Jahn- Teller systems. 
As a starting point of the model we introduce the effec- 
tive electron-configurational Hamiltonian as follows 

K 

H s = -As z - ts x - P Q + — Q 2 - aQs z , (34) 

where S2 = l(j ° x ),-. = *( J J ) are Pauli matri- 
ces, and the first term describes the bare energy splitting 



of "up" and "down" states with energy gap A (CT en- 
ergy) , while the second term describes the coupling (mix- 
ing) between "up" and "down" pseudospin states. In 
terms of a pseudospin analogy the both parameters may 
be associated with effective fields. The third and fourth 
terms in (1221) describe the linear and quadratic contri- 
butions to the configurational energy. Here, the linear 
term formally corresponds to the energy of an external 
configurational strain described by an effective strain pa- 
rameter p, while quadratic term with " elastic" constant 
K is associated with a conventional harmonic approx- 
imation for configurational energy. The last term de- 
scribes the electron-configurational (vibronic) interaction 
with a being a electron-configurational coupling constant. 
Hereafter we make use of dimensionless configurational 
variable Q therefore all of the model parameters are as- 
signed the energy units. Our model Hamiltonian has the 
most general form except the simplified form of electron- 
configurational coupling where we omit the term cx Qs x . 
In frame of adiabatic approximation the eigenvectors for 
the Hamiltonian can be written as follows: 

= cosa| t) + sina| |); 
*_(Q) = sina| T) -cosa| |), (35) 

where 

t 

tan 2a — — — . 

A + aQ 

The corresponding eigenvalues 

E±(Q) = ^Q 2 -pQ±\ [(A + aQf + t 2 ] 1/2 (36) 

define the upper and lower branches of the configura- 
tional, or adiabatic potential (AP), respectively. These 
potential curves describe the energy of |±) states as func- 
tions of configurational coordinate Q. The impact of 
different model parameters on the shape of AP can be 
easily understood, if we neglect transfer integral (t = 0) 
and deal with two ideal parabolas describing the con- 
figurational energy for "up" and "down" states | Tji); 
respectively. These identical parabolas with minima at 
QjP = (p =F § )/ K are shifted relative to each other. The 
shift in configurational coordinate is AQ = — = 
a/2K, while the shift in energy is AE^ = E { + ] - E^ = 

pa/K + A, where E^ = E ± (Q ( ± y ) or the energy of the 
bottoms of respective parabolas. Interestingly, the en- 
ergy shift is determined both by "mechanic" and "elec- 
tronic" forces. It is worth noting that the Q-shift of the 
"center of gravity" of AP is determined by the effective 
strain p. The electron-configurational coupling leads to 
a pulling apart the parabolas for "up" and "down" elec- 
tronic states. The condition 

\AE ( V\ = \pa/K + A\ =a 2 /2K (37) 
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FIG. 3: Typical adiabatic potentials of CT center: a) A = 
1.0; t = 0.1;p = -0.65; a = 1.0; K = 2.0; b)A = 0.0; t = 
0.1; p = -0.65; a = 1.0; K = 2.0; c)A = 0.0; t = 0.1; p = 
-0.65; a = 2.5; K = 2.0; d)A = 0.0; t = 0.3; p = -0.65; a = 
2.5; K — 2.0. Arrows mark CT transitions (see text for de- 
tails). 

defines two specific points, where the minimum of one 
of the branches crosses another branch, thus specifying 
the parameters range admissive of a bistability effect. If 
> E-(Q < f > ) we arrive at rather conventional situ- 
ation which is typical for long-lived CT states capable 
of decaying due to spontaneous Franck-Condon radia- 
tive transitions. In opposite case < E_(Q^?') we 
deal with only radiativclcss relaxation channels, and the 
situation strongly depends on the magnitude of trans- 
fer integral t ^ which leads to a crucial rearrange- 
ment of AP near the crossection point of two branches 
with a number of novel effects of principal importance. 
First, we obtain two isolated branches of AP, the upper 
and lower ones, respectively, with quantum superposi- 
tions ty± which describe unconventional states with a 
mixed valence of CuOm plaquettes. An illustrative ex- 
ample of typical adiabatic potentials is shown in Fig. 3. 
Upper branch of AP has a single minimum point S with 
approximately "fifty-fifty" composition | f, J.) of the re- 
spective superposition state 4 f +. 

For lower branch of AP we have either a single mini- 
mum point or the two- well structure with two local min- 
imum points (see Fig. 3c), leading to a "bistability" effect 
which is of primary importance for our analysis. Indeed, 
these two points may be associated with two (meta) stable 
charge states with and without CT, respectively, which 
form two candidates to struggle for a ground state. 

It is easy to see, that for large values of the transfer 
integral the system does not manifest bistability, i.e. the 
lower branch of conformational potential has only a single 
minimum (see, e.g., Figs. 3c, d). 

Thus we conclude that all the systems such as copper 
oxides may be divided to two classes: CT stable systems 
with the only lower AP branch minimum for a certain 
charge configuration, and bistable, or CT unstable sys- 
tems with two lower AP branch minima for two local 



charge configurations one of which is associated with self- 
trapped CT excitons resulting from self-consistent charge 
transfer and electron-lattice relaxation.^ Moreover, in the 
latter case wc deal with an additional metastable state on 
the upper AP branch which is characterized by a mixed 
valence and ionicity and seems to be intermediate one. 

Our simple model analysis allows to uncover the rel- 
ative role of different parameters governing the charge 
stability in the system. At the same time one should 
note that actually we deal with more complex multimodc 
problem with a large body of intermediate metastable 
states which can be seen in luminescence spectra. 



D. Nucleation of EH droplets and phase separation 
effects in CT unstable A/ phase 

The AP bistability in CT unstable insulators points 
to tempting perspectives of their evolution under either 
external impact. Metastable CT excitons in the CT un- 
stable M° phase being the disproportionation quanta 
present candidate " relaxed excited states" to struggle for 
stability with ground state and the natural nucleation 
centers for electron-hole liquid phase. What way the CT 
unstable M° phase can be transformed into novel phase? 
It seems likely that such a phase transition could be real- 
ized due to a mechanism familiar to semiconductors with 
filled bands such as Ge and Si where given certain con- 
ditions one observes a formation of metallic EH-liquid as 
a result of the exciton decay^ However, the system of 
strongly correlated electron M ~~ and hole M + centers ap- 
pears to be equivalent to an electron-hole Bose-liquid in 
contrast with the electron-hole Fermi-liquid in conven- 
tional semiconductors. However, the Wannier excitons 
in the latter wide-band systems dissociate easily produc- 
ing two-component electron-hole gas or plasma^ while 
small CT excitons both free and self-trapped are likely to 
be stable with regard the el-h-dissociation. At the same 
time, the two-center CT excitons have a very large fluctu- 
ating electrical dipole moment \d\ ~ 2eRuM and can be 
involved into attractive electrostatic dipole-dipole inter- 
action. Namely this is believed to be important incentive 
to the proliferation of excitons and its clusterization. The 
CT excitons are proved to attract and form molecules 
called biexcitons, and more complex clusters where the 
individuality of the separate exciton is likely to be lost. 
Moreover, one may assume that like the semiconductors 
with indirect band gap structure, it is energetically favor- 
able for the system to separate into a low density exciton 
phase coexisting with the microregions of a high density 
two-component phase composed of electron M ~ and hole 
M + centers, or EH droplets. Indeed, the excitons may be 
considered to be well defined entities only at small con- 
tent, whereas at large densities their coupling is screened 
and their overlap becomes so considerable that they loose 
individuality and we come to the system of electron M~ 
and hole Al + centers, which form a electron-hole Bosc liq- 
uid. An increase of injected excitons in this case merely 
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increases the size of the EH droplets, without changing 
the free exciton density. 

Homogeneous nucleation implies the spontaneous for- 
mation of EH droplets due to the thermodynamic fluc- 
tuations in exciton gas. Generally speaking, such a state 
with a nonzero volume fraction of EH droplets and the 
spontaneous breaking of translational symmetry can be 
stable in nominally pure insulating crystal. However, the 
level of intrinsic non-stoihiometry in 3d oxides is signifi- 
cant (one charged defect every 100-1000 molecular units 
is common). The charged defect produces random elec- 
tric field, which can be very large (up to 10 s Vcm -1 ) thus 
promoting the condensation of CT excitons and the inho- 
mogeneous nucleation of EH droplets. Deviation from the 
neutrality implies the existence of additional electron, or 
hole centers that can be the natural centers for the inho- 
mogeneous nucleation of the EH droplets. Such droplets 
are believed to provide the more effective screening of 
the electrostatic repulsion for additional electron/hole 
centers, than the parent insulating phase. As a result, 
the electron/hole injection to the insulating M° phase 
due to a nonisovalent substitution as in I^-xSr^CuO,^, 
Nd2- a; Ce a; Cu04, or change in oxygen stoihiometry as in 
YBa2Cu30e+a;, La2Cu04_5, La2Cui_ a; Li 2 ;04, or field- 
effect is believed to shift the phase equilibrium from the 
insulating state to the unconventional electron-hole Bose 
liquid, or in other words induce the insulator-to-EHBL 
phase transition. 

The optimal way to the nucleation of EH droplets in 
parent system like La2Cu04, YBa2Cu30g is to create 
charge inhomogeneity by nonisovalent chemical substitu- 
tion in Cu02 planes or in "out-of-plane stuff" , including 
interstitial atoms and vacancies. This process results in 
an increase of the energy of the parent phase and creates 
proper conditions for its competing with others phases 
capable to provide an effective screening of the charge 
inhomogeneity potential. The strongly degenerate sys- 
tem of electron and hole centers in EH droplet is one 
of the most preferable ones for this purpose. At the 
beginning (nucleation regime) an EH droplet nucleates 
as a nanoscopic cluster composed of several number of 
neighboring electron and hole centers pinned by disorder 
potential. Charged defects supporting the EH droplet 
nucleation promote the formation of metastable (" super- 
heated" ) clusters of parent phase. It is clear that such a 
situation does not exclude the self-doping with the forma- 
tion of a self-organized collective charge-inhomogeneous 
state in systems which are near the charge instability. 

EH droplets can manifest itself remarkably in various 
properties of the 3d oxides even at small volume fraction, 
or in a "pseudo-impurity regime" . Insulators in a PS 
regime should be considered as phase inhomogeneous sys- 
tems with, in general, thermo-activated mobility of the 
inter-phase boundaries. On the one hand, main features 
of this "pseudo-impurity regime" would be determined 
by the partial intrinsic contributions of the appropriate 
phase components with possible limitations imposed by 
the finite size effects. On the other hand, the real prop- 



erties will be determined by the peculiar geometrical fac- 
tors such as a volume fraction, average size of droplets 
and its dispersion, the shape and possible texture of the 
droplets, the geometrical relaxation rates. These factors 
are tightly coupled, especially near phase transitions for 
either phase (long range antiferromagnetic ordering for 
the parent phase, the charge ordering and other phase 
transformations for the EH droplets) accompanied by the 
variation in a relative volume fraction. Numerous ex- 
amples of the unconventional behavior of the 3d oxides 
in the pseudo-impurity regime could be easily explained 
with taking into account the inter-phase boundary effects 
(coercitivity, the mobility threshold, non-ohmic conduc- 
tivity, oscillations, relaxation etc.) and corresponding 
characteristic quantities. Under increasing doping the 
"pseudo-impurity regime" with a relatively small volume 
fraction of EH droplets (nanoscopic phase separation) 
can gradually transform into a macro- (chemical) "phase- 
separation regime" with a sizeable volume fraction of EH 
droplets, and finally to a new EH liquid phase. Phase 
separation is now widely discussed as an important phe- 
nomenon accompanying the high-T c superconductivity. 

Our scenario can readily explain photo-induced effects 
in pseudo-impurity phased Indeed, the illumination of 
a material with light leads to the generation of eh-pairs 
that will proliferate and grow up to be a novel nonequi- 
librium EH droplet or simply to be trapped in either EH 
droplet with the rise in its volume fraction. The excita- 
tion energy appears to be lower when exciton is closer to 
the EH droplet. Therefore once the excitation transfer 
is finite, the optical excitation is attracted to the nearest 
neighbour of the EH clusters so that this cluster expands 
effectively under the light irradiation. In other words, 
the photoexcitation would result in an increase of the 
EH droplet volume fraction, that is why its effect in op- 
tical response resembles that of chemical doping. After 
switching off the light the droplet phase would relax to 
the thermodynamically stable state. Furthermore, such 
a simple model can immediately explain the persistent 
photoconductivity (PPC) phenomena, found in insulat- 
ing YBaCuO system^ where the oxygen reodering pro- 
vides the mechanism of a long-time stability for the EH 
droplets. In PPC phenomena, an illumination of a mate- 
rial with light leads to a long-lived photoconductive state. 
During the illumination of underdoped YBaCuO near the 
insulator- metal transition, the material may even become 
superconducting. 



VI. ELECTRON-HOLE BOSE LIQUID 

Hereafter we would like to address the bosonic nature 
and some properties of the mixed-valence binary (dis- 
proportionated) M ± -phase termed as electron-hole Bose 
liquid (EHBL) in more details. For such a system, the 
pseudo-spin Hamiltonian \Tll can be mapped onto the 
Hamiltonian of hard-core Bose gas on a lattice (Bose- 
Hubbard model) which standard form can be written out 
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as follawsiW2& 

H BG =-Y^t ij (BlB j +B}B i )+J2v ij N i N j - f iJ2N i , 

i>j i>j i 

(38) 

where N, = B\B h fj, is a chemical potential derived from 

N 

the condition of fixed full number of bosons Ni = ^ (Ni), 

or concentration n = Ni/N € [0, 1]. The Uj denotes an 
effective transfer integral, Vij is an intersite interaction 
between the bosons. Here B'(B) are the Pauli creation 

(annihilation) operators which are Bose-like commuting 
for different sites [Bi, Bj] = 0, for i ^ j, and for the same 

site Bf = (Bj) 2 = 0, [Bi,Bt] = 1 - 2N t , N t = bJ-Bj; N 
is a full number of sites. 

The disproportionated phase as well as the lattice 
hard-core Bose-gas is equivalent to a system of spins 
s = 1/2 exposed to an external magnetic field in the z- 
direction. Indeed, the charge (e, h), or M^-doublet, that 
is two different valence states of M-centers, one might 
associate with two possible states of the charge pseudo- 
spin (isospin) s = i: | + |) and | — |) for electron M~ 
and hole M + centers, respectively. Then the effective 
Hamiltonian can be written as followsS^LiS 

H BG = £ JZ?(8+8j + SpT) + £ J^SIS) - hJ2*!, 
i>j i>j i 

(39) 

where J*? = 2ty, J£ = V iS , h = n - J2 Vij, s = 

-LBs+ = -^ 2 B\s* = -%+B\Bi, s ± = T^{s x ±*s v )- 
Below we make use of conventional two-sublattice ap- 
proach. For the description of the pseudospin ordering 
to be more physically clear one may introduce two vec- 
tors, the ferromagnetic and antiferromagnetic ones: 

m = l« Sl ) + (s 2 )); 1 = l«si) - ( S2 »; m 2 + 1 2 = 1 , 

where m • 1 = 0. The hard-core boson system in a two- 
sublattice approximation is described by two diagonal 
order parameters l z , m z and two complex off-diagonal 
order parameters: m± = ^-^={rn x ± im y ) and l± = 

^f-^=(l x ± ily). The complex superfluid order parameter 
\p(r) = |^(r)| exp —i<j> is determined by the in-plane com- 
ponents of ferromagnetic vector: 'J'(r) = \({Bi +62)) = 
sm_ = smj_ exp — i<j>, m± being the length of the in-plane 
component of ferromagnetic vector. So, for a local con- 
densate density we get n s = s m\. It is of interest to 
note that in fact all the conventional uniform T — states 
with nonzero ^(r) imply simultaneous long-range order 
both for modulus |^(r)| and phase <f>. The in-plane com- 
ponents of antiferromagnetic vector l± describe a stag- 
gered off-diagonal order. 

The model exhibits many fascinating quantum phases 
and phase transitions. Early investigations predict at 
T = charge order (CO), Bose superfluid (BS) and 



mixed (BS+CO) supersolid uniform phases with an 
Ising-type melting transition (CO-NO) and Kosterlitz- 
Thouless-type (BS-NO) phase transitions to a non- 
ordered normal fluid (NO) in 2D systemsi 36 i 37 i 38 



VII. TOPOLOGICAL CHARGE 
FLUCTUATIONS IN MODEL MIXED- VALENCE 
SYSTEM 

Above we focused on the homogeneous phase states 
of the mixed-valence system. Main short-length scale 
charge fluctuations in M° and systems are associ- 
ated with a thermal exciton creation, or annihilation due 
to a reaction: (M° - Af°) <-> (M+ - M~). Amongst 
the long-length scale charge fluctuations in a model sys- 
tem we would like to address the topological defects in 
quasi-2D systems such as cuprates, in particular, differ- 
ent bubble- like entities like skyrmions, or another out-of- 
plane vortices. Namely these one can play the main role 
in a nucleation of unconventional charge phases. 



A. Topological defects in A/°-phase 

The most interesting situation concerns the charge 
fluctuations in the conceptually simple monovalent M°- 
phase which seems to be a representative of traditional 
insulating oxides. From the viewpoint of the a, b- vector 
field formalism the M° system resembles, in a sense, the 
easy-axis (anti)ferromagnet. Indeed, the phase is speci- 
fied by a simple uniform arrangement of a and b vectors 
parallel to z-axis: a || b || O z , and its energy does not 
depend on the sense of these vectors. The analogy al- 
lows us to make use of simple physically clear pictures of 
(anti)ferromagnetic domain structures. 

Below we'll address two types of domain walls in M°- 
phase. The first would illustrate the long-scale fluctua- 
tion which conserves the mean on-site valence: in other 
words, (S z ) =0. The second would provide the example 
of a domain wall with a checkerboard charge ordering in 
its center. 

First of all, let note that instead of two parallel vectors 
a and b given the normalization condition, the M°-phase 
can be described by a unified vector n: a = an, b = /3n, 
and a + i/3 = exp(j«:), k € R. Hereafter, we denote 
n = njsin^cos^, sin 6* sin 0, cos 9}. Moreover, we may 
introduce a multitude of phases which differ only by the 
orientation of the unit vector n, or n-phases. Such phases 
may be considered to be the solutions of a purely bi- 
quadratic pseudo-spin Hamiltonian 3 ^ 

3 

H bq = -J2 ^({SkSiM^}^), (40) 

which can be obtained from (|19f) given the only nonzero 
"exchange" parameters Jkk with k — 4 -j-8. The quantum 
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Hamiltonian can be mapped onto classical Hamiltoniar. or for the winding number q = 1 



H bq = 2J 2 \n\ 2 J 



d 2 r 



(41) 



The constant zero value of the mean on-site valence 
(S z ) = is the common property of the n-phases. More- 
over, the mean value of all the pseudo-spin components 
turns into zero: (S) = 0. Let address the ferromagnetic 
n-phase with 180° domain walls which separate two ac- 
tually equivalent M°-domains with opposite direction of 
n vectors. The picture differs from that of conventional 
ferromagnets where 180° domain wall separates the do- 
mains with the opposite direction of magnetization. In 
the center of such a BS-type domain wall we may deal 
with a superfluid BS phase nlO z . 

In 2D systems such as cuprates it appears possible 
to form skyrmion-like topological defects like bubbles. 30 
The skyrmion spin texture looks like a bubble domain in 
ferromagnet and consists of a vortex-like arrangement of 
the in-plane components of spin with the z-component 
reversed in the centre of the skyrmion and gradually in- 
creasing to match the homogeneous background at in- 
finity. The spin distribution within a classical Belavin- 
Polyakov (BP) skyrmion is given as follows^ 



= q(p + (p : 



COS V — 



r 2<? _ y2q 

r 2 " + X 2 " ' 



(42) 



2rA 



A 2 
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2rA 

2 Ta 2 
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■ sm < 



- r 2 + A 2 ' 

(43) 

Skyrmions are characterized by the magnitude and sign 
of its topological charge, by its size (radius) A, and by the 
global orientation of the spin, or U(l) order parameter 
ifQ. The scale invariance of skyrmionic solution reflects 
in that its energy is proportional to winding number and 
does not depend on radius, and global phase (po. An 
interesting example of topological inhomogeneity is pro- 
vided by a multi-center BP skyrmionSi which energy does 
not depend on the position of the centers. The latter are 
believed to be addressed as an additional degree of free- 
dom, or positional order parameter. 



The classical Hamiltonian (|41|) has skyrmionic solu- 
tions, but instead of the spin distribution in conven- 
tional BP skyrmion 3 ^ we deal with a zero spin, but a 
non-zero distribution of five spin-quadrupolc moments 
^(4,5,6,7,8)^ or ({SiSj}) which in turn are determined 
by the skyrmionic texture of the n vector: 



((S* 2 + Sf» = 1 + cos 2 9 ; (S 2 ) = sin 2 9 ; ((S 2 - S 2 )) = - sin 2 9 cos 2(qy + <p ) ; 



(44) 



({S x S y }) =-sm 2 9sm2(q<p + cp );({S x S z }) = - sin20 cos(^ + <p ); ({S V S Z }) 

I 



sin 29 sm(q(p + </?o) 



Thus we arrive at the ring-shaped distribution of 
the effective ionicity (S 2 ) with the maximal value of 1 
(M ± -phase) at r = A, and the minimal value of at 
the ring center r = and far outside r — > oo (M°- 
phase). The "bosonic" off-diagonal complex order pa- 
rameter (V± 2 ) oc (S±) ex sin 2 9 exp(±2(qip + ip )) has a 
similar r-dependence, while the "fermionic" off-diagonal 
complex order parameter (V^) oc sin20exp(±(q(^ + < ( 5o)) 
turns into zero both at the ring center r = 0, far outside 
r — > oo (M°-phase), and at r = A, or everywhere, where 
the ionicity has a strictly definite value. 

Despite these skyrmions are derived from the toy 
model, they yield very instructive information as regards 
the probable spin texture of real solitons with BS-type 
domain walls and " superconductive" fluctuations in M°- 
phase. 

The CO -type domain walls with a nonzero mean on- 
site valence given the total J2 i (S z ) — can be ob- 
tained, in common, in a two-sublattice approximation 
with non-collinear a and b vectors. Let assume the M°- 



phase to be divided onto two sublattices, A and B, with 
&a = bA = blb = and 1^ = 1^=0. Then the 
CO-type domain wall may be represented by a gradual 
spatial non-equivalent rotation of a and b vectors in A 
and B sublattices providing the nonzero magnitude of l z 
components given I z a = ~I z b with its maximum in the 
center of the domain wall. 



B. Topological phase separation in 2D EH Bose 
liquid away from half-filling 

One of the fundamental hot debated problems in 
bosonic physics concerns the evolution of the charge or- 
dered (CO) ground state of 2D hard-core BH model 
(hc-BH) with a doping away from half-filling. Numer- 
ous model studies steadily confirmed the emergence of 
"supersolid" phases with simultaneous diagonal and off- 
diagonal long range order while Penrose and Onsager 40 
were the first showing as early as 1956 that supersolid 
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phases do not occur. The most recent quantum Monte- 
Carlo (QMC) simulations 41 * 4 ^ 4 ^ found two significant 
features of the 2D Bose-Hubbard model with a screened 
Coulomb repulsion: the absence of supersolid phase at 
half-filling, and a growing tendency to phase separation 
(CO+BS) upon doping away from half-filling. Moreover, 
Batrouni and Scalettar— studied quantum phase transi- 
tions in the ground state of the 2D hard-core boson Hub- 
bard Hamiltonian and have shown numerically that, con- 
trary to the generally held belief, the most commonly dis- 
cussed "checkerboard" supersolid is thermodynamically 
unstable and phase separates into solid and superfluid 
phases. The physics of the CO+BS phase separation in 
Bose-Hubbard model is associated with a rapid increase 
of the energy of a homogeneous CO state with doping 
away from half-filling due to a large " pseudo-spin- flip" 
energy cost. Hence, it appears to be energetically more 
favorable to "extract" extra bosons (holes) from the CO 
state and arrange them into finite clusters with a rela- 
tively small number of particles. Such a droplet scenario 
is believed to minimize the long-range Coulomb repul- 
sion. 

Magnetic analogy allows us to make unambiguous pre- 
dictions as regards the doping of BH system away from 
half-filling. Indeed, the boson/hole doping of checker- 
board CO phase corresponds to the magnetization of an 
antiferromagnet in z-direction. In the uniform easy-axis 
Z 2 -phase of anisotropic antiferromagnet the local spin- 
flip energy cost is rather big. In other words, the en- 
ergy cost for boson/hole doping into checkerboard CO 
phase appears to be big one due to a large contribution 
of boson-boson repulsion. However, the magnetization 
of the anisotropic antiferromagnet in an easy axis direc- 
tion may proceed as a first order phase transition with 
a "topological phase separation" due to the existence of 
antiphase domains. The antiphase domain walls provide 
the natural nucleation centers for a spin-flop phase hav- 
ing enhanced magnetic susceptibility as compared with 
small if any longitudinal susceptibility thus providing 
the advantage of the field energy. Namely domain walls 
would specify the inhomogeneous magnetization pattern 
for such an anisotropic easy-axis antiferromagnet in rela- 
tively weak external magnetic field. As concerns the do- 
main type in quasi-two-dimensional antiferromagnet one 
should emphasize the specific role played by the cylin- 
drical, or bubble domains which have finite energy and 
size. These topological solitons have the vortex-like in- 
plane spin structure and resemble classical, or Belavin- 
Polyakov skyrmions. 39 Although some questions were not 
completely clarified and remain open until now, the clas- 
sical and quantum skyrmion-like topological defects are 
believed to be a genuine element of essential physics both 
of ferro- and antiferromagnetic 2D easy- axis systems. 
The magnetic analogy seems to be a little bit naive, how- 
ever, it catches the essential physics of doping the hc-BH 
system. As regards the checkerboard CO phase of such 
a system, namely a finite energy skyrmion-like bubble 
domaini 4 ^ seems to be the most preferable candidate 



for the domain with antiphase domain wall providing the 
natural reservoir for extra bosons. The skyrmion spin 
texture consists of a vortex-like arrangement of the in- 
plane components of ferromagnetic m vector with the l z - 
component reversed in the centre of the skyrmion and 
gradually increasing to match the homogeneous back- 
ground at infinity. 

Recently it was shown that the doping, or deviation 
from half-filling in 2D EH Bose liquid is accompanied by 
the formation of multi-center topological defect such as 
charge order (CO) bubble domain(s) with Bose superfluid 
(BS) and extra bosons both localized in domain wall(s), 
or a topological CO+BS phase separation, rather than 
an uniform mixed CO+BS supersolid phase. 

The most probable possibility is that every bubble ac- 
cumulates one or two particles. Then, the number of such 
entities in a multigranular texture nucleated with doping 
has to nearly linearly depend on the doping. Generally 
speaking, each individual bubble may be characterized by 
its position, nanoscale size, and the orientation of U(l) 
degree of freedom. In contrast with the uniform states 
the phase of the superfluid order parameter for a bub- 
ble is assumed to be unordered. In the long-wavelength 
limit the off-diagonal ordering can be described by an ef- 
fective Hamiltonian in terms of U(l) (phase) degree of 
freedom associated with each bubble. Such a Hamilto- 
nian contains a repulsive, long-range Coulomb part and 
a short-range contribution related to the phase degree of 
freedom. The latter term can be written out in the stan- 
dard for the XY model form of a so-called Josephson 
coupling 44 i 46 i 47 

Hj = - Jjj cos(</?j -(fj), (45) 

where ipi , <pj are global phases for micrograins centered 
at points respectively, Jij Josephson coupling pa- 
rameter. Namely the Josephson coupling gives rise to 
the long-range ordering of the phase of the superfluid 
order parameter in such a multi-center texture. Such a 
Hamiltonian represents a starting point for the analysis of 
disordered superconductors, granular superconductivity, 
insulator-superconductor transition with array of 

superconducting islands with phases ipi , ipj . To account 
for Coulomb interaction and allow for quantum correc- 
tions we should introduce into effective Hamiltonian the 
charging energy^ 7 - 

H C h = -^q 2 J2 Ui ( c ^v n j ' 

hj 

where rij is a number operator for particles bound in i- 
th micrograin; it is quantum-mechanically conjugated to 
ip: m = —id/dtfi, {C~ l )ij stands for the capacitance 
matrix, q for a particle charge. 

Such a system appears to reveal a tremendously rich 
quantum-critical struct ure^S* 4 ^ In the absence of disor- 
der, the T = phase diagram of the multi-bubble system 
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implies either triangular, or square crystalline arrange- 
ments with possible melting transition to a liquid. It 
should be noted that analogy with charged 2D Coulomb 
gas implies the Wigner crystallization of multi-bubble 
system with Wigner crystal (WC) to Wigner liquid melt- 
ing transition, respectively. Naturally, that the addi- 
tional degrees of freedom for a bubble provide a richer 
physics of such lattices. For a system to be an insulator, 
disorder is required, which pins the multigranular system 
and also causes the crystalline order to have a finite cor- 
relation length. Traditional approach to a Wigner crys- 
tallization implies the formation of a WC for densities 
lower than a critical density, when the Coulomb energy 
dominates over the kinetic energy. The effect of quantum 
fluctuations leads to a (quantum) melting of the solid at 
high densities, or at a critical lattice spacing. The crit- 
ical properties of a two-dimensional lattice without any 
internal degree of freedom are successfully described ap- 
plying the BKT (Berezinsky-Kosterlitz-Thowless) theory 
to dislocations and disclinations of the lattice, and pro- 
ceeds in two steps. The first implies the transition to a 
liquid-crystal phase with a short-range translational or- 
der, the second does the transition to isotropic liquid. 
For such a system provided the bubble positions fixed 
at all temperatures, the long-wave-length physics would 
be described by an (anti)ferromagnetic XY model with 
expectable BKT transition and gapless XY spin-wave 
mode. 

The low temperature physics in a multi-bubble sys- 
tem is governed by an interplay of two BKT transitions, 
for the U(l) phase and positional degrees of freedom, 
respectively. 49 Dislocations lead to a mismatch in the 
U(l) degree of freedom, which makes the dislocations 
bind fractional vortices and lead to a coupling of transla- 
tional and phase excitations. Both BKT temperatures ei- 
ther coincide (square lattice) or the melting one is higher 
(triangular lattice) 4^ 

Quantum fluctuations can substantially affect these re- 
sults. Quantum melting can destroy U(l) order at suf- 
ficiently low densities where the Josephson coupling be- 
comes exponentially small. Similar situation is expected 
to take place in the vicinity of structural transitions in 
a multigranular crystal. With increasing the microgram 
density the quantum effects result in a significant lower- 
ing of the melting temperature as compared with classical 
square-root dependence. The resulting melting temper- 
ature can reveal an oscilating behavior as a function of 
particle density with zeros at the critical (magic) densi- 
ties associated with structural phase transitions. 

In terms of our model, the positional order corresponds 
to an incommensurate charge density wave, while the 
U(l) order does to a superconductivity. In other words, 
we arrive at a subtle interplay between two orders. The 
superconducting state evolves from a charge order with 
Tc < T m , where T m is the temperature of a melting 
transition which could be termed as a temperature of 
the opening of the insulating gap (pseudo-gap!?). 

The normal modes of a dilute multi-bubble system 



include the pseudo-spin waves propagating in-between 
the bubbles, the positional fluctuations, or quasi-phonon 
modes, which are gapless in a pure system, but gapped 
when the lattice is pinned, and, finally, fluctuations in 
the U(l) order parameter. 

The orientational fluctuations of the multi-bubble sys- 
tem are governed by the gapless XY models The rel- 
evant model description is most familiar as an effective 
theory of the Josephson junction array. An important 
feature of the model is that it displays a quantum-critical 
point. 

The low-energy collective excitations of a multi-bubble 
liquid includes an usual longitudinal acoustic phonon- 
like branch. The liquid crystal phases differ from the 
isotropic liquid in that they have massive topological ex- 
citations, i.e., the disclinations. One should note that 
the liquids do not support transverse modes, these could 
survive in a liquid state only as overdamped modes. So 
that it is reasonable to assume that solidification of the 
bubble lattice would be accompanied by a stabilization 
of transverse phonon-like modes with its sharpening be- 
low melting transition. In other words, an instability of 
transverse phonon-like modes signals the onset of melt- 
ing. The phonon-like modes in the bubble crystal have 
much in common with usual phonon modes, however, due 
to electronic nature these can hardly be detected if any 
by inelastic neutron scattering. 

A generic property of the positionally ordered bubble 
configuration is the sliding mode which is usually pinned 
by the disorder. The depinning of sliding mode(s) can be 
detected in a low-frequency and low-temperature optical 
response. 



VIII. IMPLICATIONS FOR STRONGLY 
CORRELATED OXIDES 

In this Section we suggest some speculations around 
an unconventional scenario of the essential physics of 
cuprates and manganites that implies their instability 
with regard to the self-trapping of charge transfer ex- 
citons and the formation of electron-hole Bose liquid. 



A. Cuprates 

The origin of high-T c superconductivity is presently 
still a matter of great controversy. The unconventional 
behavior of cuprates strongly differs from that of ordinary 
metals and merely resembles that of doped semiconduc- 
tor. Moreover, the history of high T c 's itself shows that 
we deal with a transformation of particularly insulating 
state in which the electron correlations govern the essen- 
tial physics. 

The copper oxides start out life as insulators in con- 
trast with BCS superconductors being conventional met- 
als. It is impossible to understand the behavior of the 
doped cuprates and, in particular, the origin of HTSC 
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unless the nature of the doped-insulating state is incor- 
porated into the theory. The problem of a doped insula- 
tor is sure much more complicated than it is implied in 
oversimplified approaches such as an effective t- J-model 
when the situation resembles that of "throwing the baby 
out of the bathwater" . 

In a case of cuprates we deal with systems which con- 
ventional ground state seems to be unstable with regard 
to the transformation into a new phase state with a va- 
riety of unusual properties. 

In our view, the essential physics of the doped cuprates, 
as well as other strongly correlated oxides, is driven 
by a self-trapping of the CT excitons, both one-center, 
and two-center liSiffi Such excitons are the result of self- 
consistent charge transfer and lattice distortion with the 
appearance of the "negative-?/" effect.^- The two-center 
excitons are associated with CT transitions between two 
CUO4 plaquettes, and may be considered as quanta of 
the disproportionation reaction 

CuO!T + CuO*T -» CuOj~ + CuOif 

with the creation of electron CuOj~ and hole CuO^ - 
centers. Thus, three types of CuC-4 centers CuO^' 6,7- 
should be considered on equal footing. In this connec- 
tion we would like to draw special attention to the lat- 
tice polarization and relaxation effects that are of pri- 
mary importance both for the formation of CT exciton 
itself and its self-trapping. It should be noted that the 
photo-excited electron-hole pairs or excitons are stabi- 
lized into self-trapped excitons (STE) accompanied with 
lattice relaxation within several pico-seconds^i 

In contrast with BaBiC>3 system where we deal with 
a spontaneous generation of self-trapped CT excitons 
in the ground state, the parent insulating cuprates are 
believed to be near excitonic instability when the self- 
trapped CT excitons form the candidate relaxed excited 
states to struggle with the conventional ground stated 
In other words, the lattice relaxed CT excited state 
should be treated on an equal footing with the ground 
state. If the interaction between STE were attractive 
and so large that the cohesive energy W\ per one STE 
exceeds the energy Er of one STE , the STE's and/or 
its clusters will be spontaneously generated everywhere 
without any optical excitation, and be condensed to form 
a new electronic state on a new lattice structured 

The minimal energy cost of the optically excited dis- 
proportionation or electron-hole formation in insulating 
cuprates is 2.0 4- 2.5 eVm^^ Interestingly, that this rela- 
tively small value of the optical gap was nevertheless used 
by Goodenough^i as argument against the "negative- 
U" disproportionation reaction 2Cu(II) = Cu(III) + 
Cu(I), or more correctly 2[CuO|~]=[CuO|~]+[CuOj _ ] 
in cuprates. However, the question arises, what is the 
energy cost for the thermal excitation of such a dispro- 
portionation? In other words, what is the energy of self- 
trapped two-center CT exciton? It is this quantity rather 
than its optical counterpart defines the activation energy 



for such a reaction. The question is of primary impor- 
tance for the self-trapping scenario. The answer implies 
first of all the knowledge of electronic and ionic polar- 
ization energies for electron and hole. The polarization 
effects with its typical energy scale of the order of sev- 
eral eV appear to be of primary importance when one 
considers different charge fluctuations in insulators. 

Several general theories of self-trapping have been pro- 
posed starting from works by Landau, Pekar and Toy- 
ozawa, however, the stages of the ST process and de- 
tailed atomistic and electronic structure of ST-excitons 
and ST-holes are still unclear even in rather simple and 
well-studied systems such as alkali-metal halides. 

As regards the STE in cuprates, we have some straight- 
forward experimental indications. A key characteristic 
of the STE is its luminescence: STE are short-lived 
luminescent states of excited crystals. The observa- 
tion of photoluminescence (PL) near 2.0 -j- 2.4 eV in 
La 2 Cu0 4 ^a near 1.3 and 2.4 eV in YBa 2 Cu 3 6! £2£^ 
near 1.78, 1.95, 2.06 eV in PrBa 2 Cu 3 0g^ is a direct 
evidence of strongly localized long-lived states related 
to self-trapped excitons or their derivatives. The near- 
infrared photoluminescence was observed in many insu- 
lating cuprates To the best of our knowledge the 
most detailed PL study was performed by Denisov et 
alml in YBa2Cu306+ K in the spectral range 0.74- 1.4 eV. 
The low-temperature PL in YBa 2 Cu306 consists of three 
peaks at 0.87, 1.07, and ~ 1.4 eV, respectively. The PL 
intensity is much stronger at small doping level. More- 
over, the doping induced PL suppression manifests itself 
more strongly for the low-energy than for the high-energy 
PL peaks. At x = 0.15 only the high-energy peak located 
at 1.28 eV (T = 10 K) survives that allows us to assume 
that the STE decay becomes more effective with doping. 
The high-energy PL peak red-shifts with the lowering the 
temperature, and its intensity decreases. 

All these features can be consistently explained in 
frames of the STE nature of PL. Different PL peak can be 
assigned both to different STE and its clusters pointing 
to the multistage character of the luminescence. 

Cuprates are believed to be unconventional systems 
which are unstable with regard to a self-trapping of the 
low-energy charge transfer excitons with a nucleation of 
electron-hole (EH) droplets being actually the system of 
coupled electron CuOj - and hole CuO® - centers hav- 
ing been glued in lattice due to a strong electron-lattice 
polarization effects. Phase transition to novel hypothet- 
ically metallic state could be realized due to a mecha- 
nism familiar to semiconductors with filled bands such 
as Ge and Si where given certain conditions one observes 
a formation of metallic EH-liquid as a result of the exci- 
ton decay^ The system of strongly correlated electron 
CuOjj - and hole CuOf - centers appears to be equiva- 
lent to an electron-hole Bose-liquid (EHBL) in contrast 
with the electron-hole Fermi-liquid in conventional semi- 
conductors. A simple model description of such a liquid 
implies a system of local singlet bosons with a charge of 
q = 2e moving in a lattice formed by hole centers. Local 
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boson in our scenario represents the electron counter- 
part of Zhang-Rice singlet, or two-electron configuration 
bf g 1 Ax g . Naturally, that conventional electron CuOj" 
center represents a relaxed state of composite system: 
"hole CuC>4~ center plus local singlet boson", while the 
"non-retarded" scenario of a novel phase is assumed to 
incorporate the unconventional states of electron CuOj~ 
center up to its orbital degeneracy. 

Thus we can introduce the concept of insulator-to- 
EHBL transition as the spontaneous condensation of self- 
trapped excitons and its clusters. However, in cuprates 
we deal with the electron/hole injection to the insu- 
lating parent phase due to a nonisovalent substitution 
as in La2- a; Sr a ;Cu04, Nd2- a ;Ce a: Cu04, or change in 
oxygen stoihiometry as in YBa2Cu306+a;, La2Cu04_a, 
L^Cui-xLi^C)^ Such a substitution provokes the nu- 
cleation of EH droplets and shifts the phase equilib- 
rium from the insulating state to the unconventional 
electron-hole Bose liquid, or, in other words, induces the 
insulator-to-EHBL phase transition. Hence, the forma- 
tion of EHBL in cuprates can be considered as the first 
order phase transition. The doping in cuprates gradually 
shifts the EHBL state away from half-filing. 

It is clear that the EHBL scenario makes the doped 
cuprates the objects of bosonic physics. There are nu- 
merous experimental evidence that support the bosonic 
scenario for doped cupratesS In this connection, we 
would like to draw attention to the little known results 
of comparative high-temperature studies of thermoelec- 
tric power and conductivity which unambiguously reveal 
the charge carriers with q = 2e, or two-electron(holc) 
transport^ The well-known relation g d ^ a — const — 
— - with \q\ = 2 1 e | is fulfilled with high accuracy 
in the limit of high temperatures (~ 700 -j- 1000-ftT) 
for different cuprates (YBa2Cu306+a;, LaaBaaCueOu-i-a;, 
(Nd2/3Cei/3)4(Ba2/3Ndi/3) 4 Cu 6 Oi 6+a; ). 



B. Manganites 

Parent manganites such as LaMnOs are antiferromag- 
netic insulators with the charge transfer gap. Fundamen- 
tal absorption band in parent manganites is formed both 
by the intracenter 02p-Mn3d transfer 59 and by the small 
intercenter charge transfer excitonSfSi which in terms of 
chemical notions represent somewhat like the dispropor- 
tionation quanta with a rather low threshold of the order 
of 3 eV, resulting in a formation of electron MnOg°~ and 
hole MnOg~ centers. The CT excitons in LaMnO*3 prone 
to a self-trapping^i and may be considered to be well 
defined entities only at small content, whereas at large 
densities their coupling is screened and their overlap be- 
comes so considerable that they loose individuality be- 
come unstable with regard to the decay (the dissociation) 
to electron and hole centers, and we come to a system of 
electrons and holes, which form an EH Bose liquid. 

An instability of parent manganite LaMnC>3 with re- 



gard to the overall disproportionation such as 

MnO^" + MnO 9 ," -► MnOg " + MnO^ (46) 

was strikingly demonstrated recently by Zhou and 
GoodenoughS The transport (thermoelectric power and 
resistivity) and magnetic (susceptibility) measurements 
showed that LaMnC>3 above the cooperative Jahn- Teller 
orbital-ordering temperature Tjt ~ 750 K transforms 
into charge-disproportionated paramagnetic phase with 
fx e ff = 5.22/^s and cooperative charge transfer of many 
heavy vibronic charge carriers. It seems rather obvious 
that with the lowering the temperature we arrive at a 
system with the well developed fluctuations of the EH 
Bose-liquid phase. Strong variation of the LaMnO*3 Ra- 
man spectra with increasing laser powes^ could be re- 
lated to the photo-induced nucleation and the volume ex- 
pansion of the EH Bose-liquid, especially, as at a rather 
big excitation wavelength A = 514.5 nm, at A = 632.8 nm 
the absorption is considerably stronger in the domains of 
novel phase than in parent lattice. 

Effective nucleation of the EH Bose-liquid in mangan- 
ites could be provoked by a non-isovalent substitution 
since this strongly polarizable or even metallic phase in 
contrast with parent insulating phase provides an effec- 
tive screening of charge inhomogeneity. Indeed, in ther- 
moelectric power (TEP) experiments with doped man- 
ganites such as Lai-^Sr^MnOs Hundley and Neumeier— 
have found that more hole-like charge carriers or alterna- 
tively fewer accessible Mn sites are present than expected 
for the value x. They suggest a charge disproportionation 
model based on the instability of Mn 3+ -Mn 3+ relative to 
that of Mn 4+ -Mn 2+ . This transformation provides excel- 
lent agreement with doping-dependent trends exhibited 
by both TEP and resistivity. 

A simplified " chemical" approach to an EH Bose-liquid 
as to a disproportionated phased naively implies an oc- 
currence of Mn 4+ and Mn 2+ ions. However, such an 
approach is very far from reality. Indeed, the electron 
MnOg and hole MnOg - centers are already the mixed 
valence centers^ as in the former the Mn valence res- 
onates between +2 and +1, and in the latter does be- 
tween +4 and +3. 

In this connection, one should note that in a sense 
disproportionation reaction (|46|l has several purely ionic 
counterparts, the two rather simple 

Mn 3+ -0 2 --Mn 3+ -► Mn 2 +-0 2 --Mn 4 +, 

and 

Mn 3 +-0 2 "-Mn 3 + -> Mn 2 +-0 1 --Mn 3+ , 

and, finally, one rather complicated 

2 --Mn 3 +-0 2 - -> 1 --Mn 1 +-0 1 -. 

Thus, the disproportionation (|46|l threshold energy has 
to be maximally close to the CT energy parameter A pc i- 
Moreover, namely this is seemingly to be one of the main 
parameters governing the nucleation of EH Bose-liquid 
in oxides. 

So far, there has been no systematic exploration of ex- 
act valence and spin state of Mn in these systems. Park et 
al& attempted to support the Mn 3+ /Mn 4+ model, based 
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on the Mn 2p x-ray photoelectron spectroscopy (XPES) 
and Ols absorption. However, the significant discrep- 
ancy between the weighted Mn 3+ /Mn 4+ spectrum and 
the experimental one for given x suggests a more com- 
plex doping effect. Subias et alS examined the valence 
state of Mn utilizing Mn if-edge x-ray absorption near 
edge spectra (XANES), however, a large discrepancy is 
found between experimental spectra given intermediate 
doping and appropriate superposition of the end mem- 
bers. 

The valence state of Mn in Ca-doped LaMnC>3 was 
studied by high-resolution Mn K(3 emission spectroscopy 
by Tyson et al& No evidence for Mn 2+ was observed 
at any x values seemingly ruling out proposals regard- 
ing the Mn 3+ disproportionation. However, this conclu- 
sion seems to be absolutely unreasonable. Indeed, elec- 
tron center MnOg 0- can be found in two configuration 
with formal Mn valence Mn 2+ and Mn 1+ (not simply 
Mn 2+ !), respectively. In its turn, the hole center MnOg~ 
can be found in two configurations with formal Mn va- 
lence Mn 4+ and Mn 3+ (not simply Mn 4+ !), respectively. 
So, within the model the Mn K(3 emission spectrum for 
Ca-doped LaMn03 has to be a superposition of appropri- 
ately weighted Mn 1+ , Mn 2+ , Mn 3+ , Mn 4+ contributions 
(not simply Mn 4+ and Mn 3+ , as one assumes in Ref . ftsjj ) . 
Unfortunately, we do not know the Mn K(3 emission spec- 
tra for the oxide compounds with Mn 1+ ions, however a 
close inspection of the Mn K(3 emission spectra for the 
series of Mn-oxide compounds with Mn valence varying 
from 2+ to 7+ (Fig. 2 of the cited paper) allows to un- 
cover a rather clear dependence on valence, and indicates 
a possibility to explain the experimental spectrum for 
Ca-doped LaMn03 as a superposition of appropriately 
weighted Mn 1+ , Mn 2+ , Mn 3+ , Mn 4+ contributions. It 
should be noted that an "arrested" Mn- valence response 
to the doping in the x < 0.3 range founded in Ref. |6l is 
also consistent with the creation of predominantly oxy- 
gen holes. 

This set of conflicting data together with a number of 
additional data 69 suggests the need for an in-depth explo- 
ration of the Mn valence problem in this perovskite sys- 
tem. However, one might say, the doped manganites are 
not only systems with mixed valence, but systems with 
indefinite valence, where we cannot, strictly speaking, 
unambiguously distinguish Mn species with either dis- 
tinct valence state. 



IX. CONCLUSIONS 

We have developed a model approach to describe 
charge fluctuations and different charge phases in 
strongly correlated 3d oxides. In frames of S = 1 pseudo- 
spin formalism different phase states of the system of 
the metal-oxide M centers with three different valent 



state M° :± are considered on the equal footing. Simple 
uniform mean-field phases include an insulating mono- 
valent M°-phase, mixed-valence binary (disproportion- 
ated) M ± -phase, and mixed- valence ternary ("under- 
disproportionated" ) M 0,± -phase. We consider two first 
phases in more details focusing on the problem of elec- 
tron/hole states and different types of excitons in M°- 
phase and formation of electron-hole Bose liquid in M- 
phase. 

Our consideration was focused mainly on a number 
of issues seemingly being of primary importance for the 
various strongly correlated oxides such as cuprates, man- 
ganites, bismuthates, and other systems with CT insta- 
bility and/or mixed valence. These includes two types 
of single particle correlated hopping and the two-particle 
hopping, CT excitons, electron-lattice polarization effects 
which are shown to be crucial for the stabilization of ei- 
ther phase, topological charge fluctuations, nucleation of 
droplets of the electron-hole Bose liquid and phase sep- 
aration effect. We emphasize an important role of self- 
trapped CT excitons in typical Mott-Hubbard insulators 
as candidate "relaxed excited states" to struggle for sta- 
bility with ground state and natural nucleation centers 
for unconventional electron-hole Bose liquid which phase 
state include the superfluid. 

Pseudo-spin formalism has appeared to be very effi- 
cient to reveal and describe different aspects of essen- 
tial physics for mixed-valence system. We show that the 
coherent states provide the optimal way both to a cor- 
rect mean-field approximation and respective continuous 
models to describe the pseudo-spin system including dif- 
ferent topological charge fluctuations, in particular, like 
domain walls or bubble domains in antifcrromagnets. All 
the insulating systems such as M°-phase may be subdi- 
vided to two classes: stable and unstable ones with regard 
to the formation of self-trapped CT excitons. The latter 
systems appear to be unstable with regard the formation 
of CT exciton clusters, or droplets of the electron-hole 
Bose liquid. The model approach suggested is believed 
to provide a conceptual framework for an in-depth un- 
derstanding of physics of strongly correlated oxides such 
as cuprates, manganites, bismuthates, and other systems 
with charge transfer excitonic instability and/or mixed 
valence. We shortly discuss an unconventional scenario 
of the essential physics of cuprates and manganites that 
implies their instability with regard to the self-trapping 
of charge transfer excitons and the formation of electron- 
hole Bose liquid. 
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